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1 Introduction 



Integrability has and continues to play an important role in the context of the AdS/CFT 
correspondence [1] , where it might provide us with the first opportunity ever to exactly solve 
an interacting field theory at the quantum level. In particular, major achievements have 
been made in the spectral problem of N = 4 SYM in the planar limit [2] by its equivalence to 
a two dimensional integrable quantum field theory; the world-sheet theory of the AdSs x S 5 
superstring in the light-cone gauge [3] . These developments conceptually culminated in the 
formulation of the thermodynamic Bethe ansatz equations [4-6] describing the spectrum 
of the world-sheet theory at finite size and finite coupling through the thermodynamics of 
an accompanying mirror model [7, 8]. 

Ultimately the structure of these thermodynamic Bethe ansatz (TBA) equations for 
the ground state stems from that of a psu(2|2) 2 invariant S-matrix [9]. Interestingly the 
centrally extended Lie algebra psu(2|2), more precisely its universal enveloping algebra, 
admits a natural quantum deformation in the sense of quantum groups [10, 11]. We will 
loosely refer to the quantum deformed algebra U g (psu(2\2)) simply as psu„(2|2). This 
algebraic structure is the starting point to construct a psu 9 (2|2) 2 invariant S-matrix, giving 
a quantum deformation of the AdSs x S 5 world-sheet S-matrix. The deformation parameter 
q can in principle be any complex number, but is typically taken to be real or a phase in 
concrete applications. In the case where q is a phase, more particularly the feth principal 
even root of unity 1 q = e*fe , the psu g (2|2) 2 invariant S-matrix [12, 13] interpolates between 
the S-matrix of the original string and mirror theory at q = 1, and in the limit of infinite 
coupling [14] the S-matrix of a relativistic semi-symmetric space sine-Gordon theory [15, 
16]. Apart from interpolating between two theories we can anticipate that the model has 
very interesting intrinsic properties, particularly when the deformation parameter is taken 
to be a root of unity. Indeed the prototypical example of a g-deformed model, the XXZ 
spin chain, shows a very interesting structure in its TBA equations when q is a root of 
unity [17, 18]; there are only finitely many Y-functions, some coupling in a non-standard 
fashion. With this in mind, in the present paper we will derive the TBA equations for the 
quantum deformed theory at q = e l fc. 

In the undeformed case of the AdSs x S 5 superstring, the derivation 2 of the ground 
state TBA equations [4-6] is based on the string hypothesis [21]. These equations imply 
the corresponding Y-system [22] with rather intricate analytic properties studied in [23, 
24]. Excited string states can be described by certain modifications of the ground state 
TBA equations [25-28], which has been explicitly done for certain string states with real 
[6, 29-31] and complex momenta [32, 33]. A particularly nice test of this approach was 
provided by the TBA equations for the Konishi operator, which were shown to agree 
numerically [34] and analytically [35] with Luscher's perturbative treatment [36-39] and 
with the dual field theory at four loops by explicit computation [40, 41], and recently at 
five loops through an approach based on four-point correlation functions of BPS operators 
[42]. Furthermore, the TBA approach has been successfully applied to the computation 

x It was argued in [16] that the S-matrix theory is most naturally defined at these points. 
2 For recent reviews in the present context see for example [19, 20]. 
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of the cusp anomalous dimension in the context of Wilson loops [43, 44], and there have 
been interesting developments in describing the spectral problem through a finite set of 
non-linear integral equations [45, 46], see also [47], complementary to the TBA equations. 
These last two topics came together nicely very recently in [48]. 

Deformations in the form of twisted boundary conditions have been previously consid- 
ered in the AdS/CFT TBA setting [49-54]. However, since twisted boundary conditions do 
not affect the mirror Bethe equations, the TBA equations are basically identical. Here on 
the other hand, the deformation affects the theory at the fundamental level of the disper- 
sion relation already, and hence we have to go through the entire procedure outlined just 
above to derive the TBA equations for the ground state of the deformed model at roots of 
unity. To go beyond the ground state we need to have a candidate asymptotic solution for 
excited states, and we will address this question in an upcoming publication [55]. 

The p5U 9 (2|2) 2 invariant S-matrix which lies at the basis of our derivation of course 
satisfies all usual requirements of integrability. However, its physical unitarity has not 
been fully investigated before and as such remained an open question. We have found 
that this S-matrix (as well as the psu,j(2|2) invariant i?-matrix of course) is not physically 
unitary but rather physically pseudo- unitary 3 , in particular in a local fashion 4 . Now in 
general pseudo-unitary (pseudo-Hermitian) models come in two classes; the Hamiltonian 
has a self-conjugate or a real spectrum. In the latter case the model is quasi-unitary 
(quasi-Hermitian) [57]. Our uniformized two body S-matrix generically satisfies generalized 
pseudo-unitary for complex arguments, and has a unitary spectrum for both the string and 
mirror theory. This quasi-unitary structure does not appear to be compatible with locality 
however, and as such does not necessarily extend to the many-body scattering theory. In 
fact, the many-body S-matrix makes a clear distinction between the string and mirror 
theory; the string theory many body S-matrix has non-unitary eigenvalues and is therefore 
only pseudo-unitary, while the mirror theory many body S-matrix appears to remain quasi- 
unitary. This difference means that we can expect the mirror theory to have a stronger 
sense of reality than the string theory, and in fact this is what we observe as we will discuss 
soon; these properties translate to the string hypothesis and TBA equations. 

In general the question of a g-deformed TBA at roots of unity is notoriously difficult 
to answer. It might be surprising to learn that the TBA for a general su q (N) spin chain 
is in fact not known. One of the difficulties is related to the fact that the associated 
local Hamiltonian is non-Hermitian 5 for N > 2, highlighting the special status of the XXZ 
model. For su q (3) work has been done on the TBA for complex Toda theories [58], involving 
interesting but rather unusual non-unitary string complexes; still the general story remains 
unknown to the authors' current knowledge. Nevertheless, likely due to the quasi-unitarity 
of the mirror theory, we appear to be able to correctly account for the thermodynamics of 

3 The notion of pseudo-unitarity is defined in appendix A, see eqn. (A. 11). For the notion of pseudo- 
Hermitian quantum mechanics, see e.g. [56]. 

4 In other words, the Hermitian automorphism involved in the pseudo-unitarity factors over the one 
particle states, cf. eqs. (A. 11) and (A. 12). 

5 By a quick investigation they (naturally) appear to be pseudo-Hermitian, but not quasi-Hermitian, 
perfectly in line with our discussion. 
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the quantum deformed Hubbard model 6 with a 'real' string hypothesis, giving an elegant 
structure highly analogous to the XXZ case. Thus it appears that the mirror su g (2|2) 
theory is the next simplest case beyond the XXZ model. The coupling between the two 
copies of the deformed Hubbard model merges with this structure very neatly. 

As we will show below, the spectrum of excitations of the deformed mirror model in 
the thermodynamic limit is severely constrained. This is analogous to what happens in the 
XXZ spin chain at roots of unity, though more involved. As a result, we obtain a finite set 
of 5k TBA equations 7 which couple in an intricate way. 

V, 



6 




Figure 1. The TBA structure of the XXX model and the XXZ model at roots of unity. The XXX 
model has a standard nearest neighbour coupling between infinitely many Y-functions (the vertices 
in the diagram), while the XXZ model has only a finite number of them. In the XXZ model the 
last two Y-functions are related directly algebraically (represented by the lime-green ellipse) and 
only couple back to the next to last regular Y-function. 



The XXZ model with k = 2 is equivalent to free fermions, and also in our model the 
point k = 2 is degenerate 8 , therefore we will consider k > 2 in the present paper 9 . With 
q = e l k and k > 2, the string hypothesis for the XXZ spin chain allows only k different 
types of string solutions, and the standard TBA equations are cut at level k — 2 while 
adding an end structure between two algebraically related Y-functions and the last regular 
one [17, 18], see figure 1. Of course the algebraic relation can be used to eliminate one of the 

6 The psu ? (2|2) invariant R-matrix of [10] gives models which are closely related to the Alcaraz-Bariev 
model [59], however the detailed relation is not yet fully understood in all cases. The Alcaraz-Bariev model 
can be viewed as a quantum deformed version of the Hubbard model and is sometimes referred to as such. 
We loosely refer to our psu 5 (2|2) invariant model as the quantum deformed Hubbard model here. 

7 3fc for left-right symmetric states. 

8 For instance, on top of the obvious degeneration of the Bethe-Yang equations, the shape of the physical 
region on the torus degenerates dramatically. 

9 There is no principal objection to treating the case where q is a generic phase along the lines of the 
XXZ model [17, 18], but already there the technicalities become somewhat involved which would lead to a 
truly undesirable level of technical details here. 



- 4 - 



last two Y-functions in exchange for a nonstandard coupling; this pictorial representation 
is chosen to emphasize their origin as two separate Y-functions. 



• • • • 6 




Figure 2. The TBA structure of the Hubbard model (left) and the deformed Hubbard model at 
roots of unity (right). The green, teal (with ±) and blue dots indicate what we call Y M i w , Y±, and 
Ym\vw functions respectively; their coupling is nearest neighbour apart from the indicated coupling 
near the corner. The modification of this structure in the deformed case is analogous to the XXZ 
case on each of the su 9 (2) wings. (The reflection of the diagram is for purely aesthetic purposes.) 



For our model we will firstly need to identify the TBA structure of the g-deformed 
Hubbard subsystem. The undeformed Hubbard TBA structure is that of two su(2) wings 
coupled via two extra Y-functions. Quite naturally, the deformation results in an XXZ-like 
modification of each of these wings independently, giving the structure illustrated in figure 
2 at roots of unity 10 . This means that for the deformed Hubbard model at roots of unity 
we have 2k + 2 Y-functions with two algebraic relations between them, resulting in 2k 
independent equations. 

Coming to the full AdSs x S 5 model then, in the undeformed case the Y-system [22] is 
given by coupling the Y-systems of two Hubbard models through physical mirror excitations 
and implies the TBA structure indicated in figure 3. 

As we will show in this paper, when we analogously couple the two psu 9 (2|2) systems 
the spectrum of momentum carrying particles of the deformed mirror theory naturally ter- 
minates at bound states of length k. These bound states and their scattering properties 
arrange themselves perfectly with the subsystems, resulting in the TBA structure illus- 
trated in figure 4. There are 5k + 4 Y-functions with four algebraic relations between four 
pairs of them, resulting in 5k independent TBA equations. 

The deformed model shows various interesting features on top of the fact that it can be 
described through a finite number of TBA equations. To start with, as already indicated 

10 The TBA diagram for the generically deformed Hubbard model would have infinite extent like the 
undeformed case. 
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Figure 3. The structure of the AdSs x S 5 TBA. Two copies of the Hubbard model are connected 
via the momentum carrying nodes (purple) of the mirror theory, corresponding to Ig-functions in 
our notation. Note that Y\ couples to F_ , but no Yq couples to Y + in a local fashion. 






Figure 4. The structure of the full deformed Y-system. The coupling between the two deformed 
Hubbard systems is via k momentum carrying YQ-functions, where the kth one precisely couples 
to the end structure of the deformed Hubbard model, notably in a local fashion. 



the S-matrix upon which the TBA equations are based is not unitary but pseudo-unitary. 
Nonetheless the TBA equations are real, which is no doubt due to the not only pseudo but 
apparent quasi-unitarity of the mirror theory, briefly discussed in appendix A. At the level 
of kinematics, the physical region of complex rapidities of the mirror theory is deformed and 
in particular the real mirror line on the rapidity torus is now a true interval, as explained 
in the text around figure 6. Next, as a cute result of the periodicity introduced by the 
deformation it is possible to define a mirror x function which maps a single complex u- 



plane onto half of the torus in one go, compared to the quarter of the torus covered in the 
undeformed case, cf. eqs. (2.17) and (2.18) and figure 6. This allows for a description 
of the crossing transformation on the u-plane directly, if desired. Furthermore, the end 
structure at the level of the auxiliary TBA equations alluded to above shows an interesting 
coupling. The last two Y- functions introduced on either of the su q {2) wings turn out to be 
exactly inverse to one another (eqs. (5.4) and (C.25)), and result in a doubled contribution 
to the last regular simplified TBA equations and associated Y-system as shown in eqs. 
(6.2,6.3,6.11) and (6.12,6.16,6.23), respectively. The fact that it is possible to obtain a 
Y-system equation for length k mirror bound states is due to quite nontrivial cancellations 
around this boundary, in particularly involving the g-deformed dressing phase as discussed 
in section 6.1. These cancellations are moreover precisely such that in the end we have a 
doubled contribution from the vw functions resulting in a boundary Y-system equation of 
the form 

v + v n O+f^-V 

Y+Y k _ Q =i,2 v r fc-n™/ 
~Y*—- 1 + YkTi ' 

The solution of the ground state TBA equations shows that the ground state of the de- 
formed 'string' theory has zero energy. 

This paper is organized as follows. In section 2 we introduce our main variables 
which parametrize the fundamental representation of the q-deformed algebra, and discuss 
the structure of the underlying rapidity torus and the dispersion relation of the g-deformed 
model and its mirror. Then in section 3 we present the Bethe-Yang equations of the mirror 
model, followed by a discussion of two- and multi-particle mirror bound states in section 
4. Based on this analysis and the analysis done in appendix B we proceed to formulate the 
string hypothesis for our model in section 5 and derive the canonical TBA equations. We 
find their simplified form and associated Y-system in section 6. Finally, in section 7 we use 
the TBA equations to compute the energy of the ground state and show that it vanishes 
in complete analogy with the undeformed case. These are the main results of our paper. 
In the conclusion we recapitulate our findings and indicate interesting open questions as 
well as some preliminary findings not addressed here. We would like to emphasize that 
our string hypothesis in principle contains the string hypothesis of the quantum deformed 
Hubbard model at roots of unity, which has not been investigated before. Technical details 
regarding the S-matrix (including its pseudo-unitarity), string hypothesis, TBA kernels, the 
simplified TBA equation for Y/-, and representation theory of psu g (2|2) have been collected 
in five appendices. 

2 Kinematics of the g-deformed model and its mirror 

Here we present various useful parametrizations of the fundamental representation of the 
centrally extended psu g (2|2) algebra as well as the dispersion relations of the g-deformed 
model and its mirror. For a detailed discussion of (atypical) representations of psu 9 (2|2) 
with q a root of unity we refer the reader to appendix F. 
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2.1 Rapidity torus, u- and n-planes 

Let us start by recalling that the basic variables x^ parametrizing a fundamental repre- 
sentation of the (/-deformed algebra satisfy the following constraint [10, 11] 

K-+^)-4-+^H)H)' 

where the parameter £ is related the coupling constant g as 

P = -j - ^-^-1 (2.2) 

2 y/l-i(q-<T 1 ) i 

As in the undeformed case the fundamental variables x can be uniformized on an elliptic 
curve. This elliptic curve has real period 2oj\{k) = 4K(m) and imaginary period 2u>2( K ) = 
4«K(1 — to) — 4K(m), where K(to) is the elliptic integral of the first kind considered as a 
function of the elliptic modulus m = k 2 . It turns out convenient to introduce a variable zq 
in place of the deformation parameter q, related to it as 

= e iam(2^ ) = CS + idll _ ( 2 .3) 

cso — idno 

Here and below we use a concise notation for Jacobi elliptic functions; no subscript denotes 
the free variable in the equation, while the zero subscript indicates evaluation at zq. For 
instance, cs(#o) = csq. The modulus k is related to the coupling constant g as 



9 = ~Jl-^ - (2-4) 

In the limit q — > 1, i.e. zq — > 0, we recover the familiar relation k 2 = —4g 2 . With these 
conventions, the variables x^ are the following meromorphic functions on the z-torus 

, 2 cs + cs o dn — icso _ 1 dn + dno cs — zdno 

x + z =Ksn - — , x {z) = 2 a : — ' ( 2 - 5 ) 

an — dno cs — ^ano k sn^ cs — csq dn — icsq 



while £ is given by 



C ■ SI1 Cn / oft N 

i = -in— . 2.6 

dn n 



In this paper we want to take q to be a root of unity and therefore restrict ourselves to 
real zq. Moreover, for q to cover the unit circle once, the corresponding values of zq can be 
restricted to the domain 

- ux (k) /2<z < ui (k) /2 . (2.7) 

We further take the modulus k to be purely imaginary with Iitik > 0, so that the coupling 
constant g is positive. Accordingly, £ lies between zero and one. Note that for real z$ 
conjugation of the x^ functions on the torus takes the form 

x ± (z)* = x T (z*). (2.8) 
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Figure 5. The division of the torus by the curves \x ± \ = 1 (blue, dashed) and lm(x ± ) = 
(two greens, solid) for five different values of q on the unit circle. The orange dot on each plot 
indicates the position of the corresponding zq; the middle plot corresponds to zq just below the 
value for which q = i. The shape of the 'string' region (lar^l = 1) smoothly changes from the "fish" 
to the "hour-glass" shape as Zq runs from to Wi(/c)/2, i.e. q from q = 1 to q = —1, but it does 
interchange the location of its crests and throughs at q = i. For the 'mirror' region (lm(x ± ) = 0) 
the change around and the situation at q — i is much more dramatic as the figure shows. 



For reasons already indicated in the introduction, we will take q = e % k with k an integer 
greater than two. The special nature of k = 2 can be seen in figure 5, which show that 
the 'string' and 'mirror' type regions 'flip their character' precisely at q = i. The central 
charges of the fundamental representation are 

V 2 _ l x+ +£ _ CS + idn _ e i( a m(z+z )+ a m{z-z )) ^ ^.g) 

q x~ + £ cs — idno 

y-2 = q * + + £ = CS + idn = e i(am( 2+2o )-am(,-,o)) ( 2 1Q) 

x~ x + + £ cso — idn 
and the shortening condition reads 



9 (1- U 2 V 2 )(V~ 2 - U- 2 ) = 1 . (2.11) 



q-q- 1 J 4 

To proceed, let us introduce a multiplicative evaluation parameter u = u(x) defined as 



x+l+C+l 



(2.12) 



so that u(x + ) = q 2 u(x~) and x^ can be determined in terms of a single variable x as 

u(x + ) = q u(x) , u(x~) = q^ 1 u(x) . (2-13) 

This evaluation parameter has the following expression via the torus variable z 

csodno + cs dn 



csodno — cs dn 



(2.14) 
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From eq.(2.12) we determine the inverse function x(u) as 



x u 



U + l\ 1 ' " ! ' 



(2.15) 



There are two branch points < j^l < 1 and 1 < j^l since £ is positive and less than 1. 
The map x(u) is chosen in such a way that the cuts are ] — oo, ^|] U [^|, +oo[. As u runs 
over the interval < u < ^| its image x(u) spans the unit half-circle in the lower half 
plane. Correspondingly, l/x(u) spans the upper half-circle. Note that the variables x^ are 
expressed via x(u) in a simple manner 11 

x ± = x{q ±l u) . (2.16) 

This x function maps the u-plane onto the green regions of the torus indicated in figure 6. 
As the branch cuts of x(u) are straight lines, the branch cuts of ar^u) necessarily intersect 
and in this fashion cut of part of the u-plane. It is precisely the disconnected green region 
on the torus which corresponds to this disconnected region of the u-plane. 

For the work that follows it will provide fruitful to work in an additive setting 12 , which 
is obtained by the u-plane parametrization u = q~ %9U = e~k~ . Since the shift u — > u + — 
leads to the same value of u, the u-plane is an infinitely-sheeted cover of the u-plane. For 
this parametrization the map x(u) turns into 

e^t ( sinh ^ - i ,/ 5 2 sin 2 f -sinh 2 ^) - g 2 sin 2 f 
x{u) = ^ 1 . >- , (2.17) 



ffsinf/L + sWf 



■Kgu Trgu 

where we have replaced which originally appears under the square root by in 
front of it, removing a square root ambiguity. Because of this the function is no longer 
periodic with period 2ik/g, but rather 

" + f) = ^ < 2 - 18 > 

meaning x(u) is periodic on the u-plane with period In other words, by resolving a 
square root ambiguity, we have extended our x function beyond the original mirror theory 
u-plane. Extended in this fashion, the mirror x function covers a full vertical band of the 
torus, as illustrated in figure 6. In fact, eq.(2.18) is nothing but the crossing transformation. 
Moreover, the scattering matrix of the g-deformed model, see appendix A, can be put on 
two copies of the u-plane, where it is compatible with crossing symmetry and its matrix 
part is periodic with period 4ik/g in either argument. 



11 The parameters x are two complex variables obeying one (complex) relation (2.1). Expressing x in 
terms of u, this relation is explicitly resolved. 

12 We will shortly see that this is exactly the analogue of the conventional hyperbolic parametrization of 
the XXZ model with |A| < 1. 
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k/i 



-k/g 



-2k/g 



" — 





Im x + > 0, Im x 


>0, 


Im x + < 0, Im x 


<o. 




Im * + > 0, Im x 


<0, 
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Im x + < 0, Im x~ 


>o, 




|* + |>1, \x~\>\, 




IB 


\x + \>\, \x~\<\, 






\x + \<l, \x'\>\, 






\x + \<l, bTkl. 





Figure 6. The torus and the mirror u-plane. The left figure shows the torus with the distinguishing 
properties of the x ± functions defined on it. The right figure depicts the mirror u-plane, with the 
same indications regarding x ± , illustrating the map between the mirror u-plane and the torus. 



The x function has branch cuts on the u-plane running outward along the real line 
from ±U{, where 

k . 1 + £ 2k ( . ir\ . 

u\j= — log = — arcsmh g sin — , (2-19) 

Trg 1 — £ Tig \ k) 

as well as outward from the points ±Ub + 2ik/g as illustrated in figure 7. In the limit 




Figure 7. The absolute value and imaginary part of x(u) on the u-plane. The green surface 
distinguishes \x\ > 1 from |x| < 1 in the left plot, while the blue surface distinguishes the sign of 
the imaginary part in the right plot. 
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k — > oo, the variable x(u) tends to the standard mirror variable of the undeformed theory 

x(u) = -(n — i\J 4 — u 2 ) , (2.20) 

and naturally has the conjugation property of the standard mirror variable; x(u)* = x ^*^ ■ 

2.2 The dispersion relations 

It is natural to make the following identification 

V = q% , U = e^ , (2.21) 

where H and p are the energy and momentum of the model. For real z, the energy is a 
non-negative periodic function, while the momentum takes values in the interval (— vr, ir) 
as z runs over (— -^p) as in the undeformed case. In the limit q — > 1 we smoothly obtain 
the undeformed string theory result p = 2amz and H = dnz. Note that in terms of H 
and p thus introduced, the shortening condition turns into the following dispersion relation 

e 2 [H/2)l-g 2 s^ P - = [l/2} 2 q , (2.22) 

where we have introduced a manifestly positive "coupling constant" e 

e 2 = l- ^(q-q- 1 ) 2 = l+g 2 S m 2 ^-, (2.23) 



and the notation 

J 9 



•1-7^' < 2 - 24 > 



To obtain the dispersion relation for the mirror theory, we perform the double Wick rota- 
tion, which just as in the undeformed case amounts to the following replacement in (2.22) 

H^ip, p^iH, (2.25) 
where p and H are the momentum and energy of the mirror model. In this way we find 

H = 2 arcsinhM^l _ e 2 [ip] 2 1/2 . (2.26) 

In particular, for q = e~k the formula takes the form 

(1 sin \ _ sinh 2 ^p \ 
— r^fi 1 + e 2 5-irH ■ 2-27 
g sin f y sin 2 f k J 

Very importantly, the mirror momentum is real on the line at ^ i n ^ ne green region of 
the torus, but not on this line in the yellow region. This is perhaps surprising, but most 
certainly not a problem since this interval on the torus already covers the whole real line 
of mirror momenta; this interval corresponds to the whole real line of the u-plane. On this 
same interval the energy is positive and bounded from below. There are two natural ways 
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to consider the limit g — > oo. We can take the limit g — > oo with k and p fixed, resulting 
in a linear (phononic) dispersion 

H=\\p\. (2.28) 



> „„>, „„„ u,, <„, . £ »^ *£ 

[13], we obtain the relativistic dispersion 



Alternatively if we first rescale the energy and momentum as H — > — and p — > - 1 as in 



H 2 -f = — W ( 2 - 29 ) 



cos 



2 jr 
2k 



up to a rescaling of g by a factor of two. Finally, let us note that the mirror transformation 
on the rapidity torus is given by a shift of the z variable by a quarter of the imaginary 
period, as in the undeformed case. 

3 Bethe-Yang equations for the mirror model 

We are interested in constructing the g-deformation of the TBA equations for strings on 
AdSs x S 5 . Thus we assume that the corresponding S-matrix has psu g (2|2) © psu g (2|2) 
symmetry and therefore factorizes into two copies each invariant under psu g (2|2). As in 
the undeformed case, this uniquely fixes the S-matrix up to a scalar factor. Of course 
this scalar factor is not uniquely determined by the requirements of unitarity and crossing, 
but a very natural generalization 13 of the undeformed scalar factor that satisfies these 
requirements has been found [12]. We will take this as the S-matrix defining our theory 
and proceed from there. 

The auxiliary problem corresponds to two copies of the quantum deformed Hubbard 
model, and analyticity of the associated transfer matrix 14 implies the following Bethe-Yang 
equations for a single copy of the auxiliary roots y m and w n 

k 1 _ r^K™ 

l = UV« yi ^\^Il YlI ^ i , (3-1) 

K 11 K™ o 

±A : w n - qvi "~ q z w n - wj 
i=i j=i J 

In the n-plane parametrization 

■kqu ngv irgw 

u = e~ , v = e~ , w = e~k~ , (3.3) 



13 Taking a diflterent scalar factor corresponding to a different solution to the crossing equation would 
change the S-matrix. However, insisting on the proper q — > 1 limit we expect the modification not to 
change the bound state picture we describe below. Therefore the only modification would be a simple 
prefactor in the s[(2) S-matrix which should not affect our results in any further way. 

14 This transfer matrix will be explicitly presented in an upcoming paper. 
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the auxiliary Bethe equations take the form 



X — TT ~ x i . x t TT ™ { — ~ Wl ~° } (34) 

1 _ S 5inh 1 K i^+j) t~t sinh I ^ ~ W A - f ) ; (3 . 5) 
'=' sinh 1 K - ^ - f ) 7=i sinh 1 K - wj + 1 ) ' 



where 



e^= v = J j 5-. (3.6) 

As mentioned above, the second set of auxiliary equations is identical to those of the 
inhomogeneous Heisenberg XXZ spin chain; the limit V{ — > gives the homogeneous XXZ 
spin chain. In the limit k — > 00 with Bethe roots kept finite, by construction we get the 
auxiliary Bethe- Yang equations of the undeformed model 



(3.7) 
(3.8) 




i=l 11 1 g j=l " 



where v = y + i. With the definition of momentum (2.21,2.25) of the previous section and 
the scalar factor of the S-matrix, we can write down the full set of Bethe- Yang equations 
of the quantum deformed mirror model 

i=^4%fe^npfei? (3 - 9) 

a=l i=l Hi x l V I 




4- (a) ■ 1_ 7TO / (a) (a) 

w — TT 7^ — — . 3 - 10 

sinh 1 (u4 a) - + j) ^ sinh f (wi a) - wf - f ) 

TT — 1 \ — — — TT — < — 7-1 — — . 3 -ii 



sinh 1 (uT - <T " |) ,ti ^nh ff (^ - < J + 



9 



where a = 1,2. Here, as usual, p/ stands for the mirror momentum of Zth particle and R is 
a length of the mirror circle. We will come back to the scalar factor S S U2) when discussing 
bound states in the next section, and have summarized the full deformed S-matrix and its 
properties in appendix A. 

4 Bound states of the mirror theory 

In order to discuss the thermodynamics of the mirror theory in the infinite volume limit, we 
need to determine the spectrum of excitations that make up the thermodynamic ensemble 
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in infinite volume. In this section we discuss the spectrum of physical excitations of the 
infinite volume mirror model, while in the next section we join these with the spectrum of 
auxiliary excitations whose analysis we have carried out in appendix B. Together this goes 
under the name of the string hypothesis. 

4.1 Two-particle bound states 

To discuss the physical bound states of the mirror theory we need to analyze consistency 
of the Bethe equations in the limit R — )■ oo. In the absence of auxiliary roots, the main 
Bethe equation (3.9) takes the form 

K l 

l = e^ R l[S m (x h x i ), (4.1) 

where the S-matrix corresponding to the ^-deformed analogue of the st(2) sector of string 
theory is given by 

l-^-f sinhff (ui - u 2 + f ) f 1 ' 



x - - x + 1- - l— sinh - !i 2 - f ) \ J ■ (4 ' 2) 



and a is the dressing factor. An explicit formula for p(u) is given below in eq. (4.13). For 
complex values of momenta the S-matrix exhibits a pole at = ■ At the level of the 
Bethe equations this pole is accompanied by a divergence of the momentum factor in the 
limit R — > oo. Indeed, if the first particle has a momentum with positive imaginary part, the 
factor e iplR in (4.1) goes to zero, but the total equation remains finite when accompanied 
by a pole at = . If the second particle has a momentum with negative imaginary 
part such that the resulting total momentum pi +p2 is real, we obtain a two-particle bound 
state. 

The two-particle bound state condition can be solved in terms of the torus z- variable 
and the solution corresponding to the simplest two-particle bound state is shown in figure 
8. As in the undeformed model, there is a critical value of the particle momenta; below the 
critical value the rapidities of the constituent particles are conjugate to each other (with 
respect to the real line of the mirror theory) while beyond the critical value the conjugation 
property is lost. Because of this, in addition to the solution indicated in figure 8 there is 
another solution corresponding to reflecting the picture in the mirror line. This behavior 
is completely analogous to the undeformed case [8], and has been observed in the present 
context in [13]. 

The critical value of the momenta for the two-particle bound state is given by p cr = 
P2(ub). Explicit computation gives 

k. cos 2 f + g 2 sin 2 f + Pl sin 2 f + p 2 sin f 

Per = ~ fog * , , 2 ■ 27T ' 43 

vr 1+5 sm z | 
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Figure 8. The two-particle bound state on the torus with its momentum and energy. The con- 
stituents of the bound state lie along the orange curves, where p = p% +p2 increases from — oo to oo 
as we move from point A to B. The kinks corresponds to critical values of the momentum, beyond 
which the two constituent momenta arc no longer complex conjugate. This is shown in the upper 
right graph where the green and yellow lines show the real (solid) and imaginary (dashed) parts of 
the constituent momenta. 



where 



in - , 1 • V • b/ 'shr J\ (4.4) 

7T TT 7T 7T TT 

p 2 = A /-2cos 2 - + 2g 2 sin 2 - + 4g 4 sin 4 - + 2(cos 2 - + g 2 sin 2 -)pi . (4.5) 

* ' Zr* k* i"* Zr* k* 

IV l\i Ft IV IV 



In the limit k — > oo we get p cr = y — 2 + 2^/l + 4g 2 which coincides with the result 
obtained in [8]. 

Note that the two-particle mirror bound state trajectory, together with part of the 
two-particle string bound state trajectory which lies along the boundary of the small blue 
region in figure 8 isolates a region of the torus which can be naturally called the physical 
mirror region of the theory. 

4.2 Multi-particle bound states 

In case the momentum of the two-particle bound state is not real we will necessarily have 
to involve a third particle to render the first two equations finite; the bound state grows. 
Without loss of generality we can take p\ + p~2 to have positive imaginary part 15 , giving a 



Otherwise we could have equivalently analyzed the poles of S starting with the Bethe-Yang equation 
for the second particle. 
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zero that is cancelled by the pole 3 . At this point either the total momentum is 

real again, giving a three particle bound state, or we add a fourth particle and continue 
the process. In this way, we can obtain a Q-particle bound state defined by the conditions 

x l = x 2 > x 2 = x 3 > • • • ' X Q-l =x Q- (4-6) 

The question now becomes which bound states are physical. 

As in the undeformed case, we can insist that the bound state condition has a unique 
solution in the physical region of our theory. The physical mirror region referred to just 
above, the big green region of figures 6 and 8, has this property, and contains bound states 
length of up to and including k. Were we to go beyond length k we would first enter the 
lower blue and upper olive-green regions of the torus, as will become clear shortly when 
we consider bound states on the tt-plane (c/. eq. (4.7) just below). These regions contain 
in particular fundamental particles and on their boundaries two particle bound states of 
the string and anti-string theory, making it undesirable to consider proceeding into them. 
Were we to continue even further, we would enter what we can by now refer to as the anti- 
mirror region and we would manifestly lose the uniqueness of the solution to the bound 
state equation. Moreover, this would correspond to constructing bound states containing 
both particles and anti-particles. 

For these reasons we define the physical region of the mirror theory as the large green 
region on the torus of figures 6 and 8. On the u-plane this corresponds to the strip 
|Im(u)| < k — 1, as is clear from figure 6. 

For a more concrete discussion, let us put our bound states on the u-plane. There 
the set of roots making up a Q-particle bound state takes the standard form of the Bethe 
string 

Uj =u + -(Q + l-2j), j = l,...,Q. (4.7) 

The pole structure of these bound states is compatible with the Bethe- Yang equations since 
the mirror momentum p{u) has positive imaginary part in particular for < Im(u) < k — 1 
as illustrated in figure 9. This means that the imaginary part of the momentum of the first 
half of the particles is positive for all bound states we consider 16 . 

Hence in the limit R — > oo we have Q-particle bound states of length up to and 
including k defined by the following equations 17 

x i = x t ' x 2 = x t ' ■ • • ' x Q-i = X Q ■ (4-8) 

16 At this point the attentive reader might be slightly concerned that by these considerations higher length 
bound states also appear to be allowed, at least as far as the imaginary part of the mirror momentum is 
concerned, while we just argued they should not be. The only reason we are even pondering this question 
at this point is that due to the deformation the it-plane can naturally cover both the mirror and anti-mirror 
region, making it possible to smoothly consider the mirror momentum along the full imaginary direction of 
the torus in one go. In the undeformed limit, we could construct a very similar picture, but (un) fortunately 
it can only be obtained by gluing two entire planes together. Colloquially speaking, fusing through the line 
ik/g is like fusing through infinity in the undeformed case. 

17 Note that due to the periodicity introduced by the deformation, this type of string configuration for 
Q = k necessarily implies xf = 1 /x^ . This appears to contradict the string hypothesis since the improved 
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Im(u) 



Figure 9. The left plot shows the imaginary part of p(u) in blue for some representative value of 
Re(u), the thin green line shows \e^\. The shaded area of the plot lies outside the physical mirror 
region. The right plot shows the same functions for the undeformed model. The maximal value of 
Imp in reached at Im(u) = k/g in the deformed model and at Im(u) = oo in the undeformed one. 



These Q-particle bound states transform in short representations of the symmetry algebra 
with central charges 



U 



i x+ + i 

qQ x~ + £ 



V6 = <t 



X x + + £ 



Here the variables x^ = x(u ± iQ/g) satisfy the relation 



1 V+ 1 

x~ 



x + 



1 



£ + 



The associated mirror momentum and energy are given by 

p Q = -£io g ^, i Q = -]o g u% 



(4.9) 



(4.10) 



(4.11) 



Both the energy and momentum are real quantities for real values of u, as follows from the 
conjugation property of x(u) 



x + {u) 



On the u-plane they are explicitly given by 



1 



x (u*) 



(4.12) 



PQ 



log 

77 coshff(u 



cosh + — i\l q 2 sin 2 £ — sinh 2 ^ 



S Q = - log 



sinhf^ + f) 



5 2 sin 2 I- sinh 2 =f(«-f 



2 sin 2 f-sinh 2 l(n+f 



sinh l(^- f 



(4.13) 
(4.14) 



1 \ I X< X 

dressing phase Yl(xi,Xk) = _ % cr(%i,%2) has an apparent zero there. However, as the name suggests 

the improved dressing phase is perfectly finite at these points; the zero is precisely cancelled by a pole in 
a(xi, x 2 ). 



-18- 



5 The string hypothesis and TBA equations 

In this section we describe the spectrum of the infinite volume mirror model and derive 
the TBA equations describing thermodynamics of the mirror theory. 

5.1 Physical excitations of the mirror model 

Based on the above analysis of bound states of the mirror theory and the analysis carried 
out in appendix B, we propose that the spectrum of excitations of the g-deformed mirror 
model in the thermodynamic limit is given by 

1. Bound states of fundamental particles; Q-particles, 

2. Complexes of y and w-roots; M|uti;-strings, 

3. Complexes of w-roots; M|w;-strings, 

4. Single y-roots; y-particles. 

The significant difference with the string hypothesis for the undeformed string mirror model 
[21] is that the length and type of such complexes is considerably constrained; we have put 
the technical discussion appendix B and present an overview here. 

Q-particles are bound states of fundamental particles, which in the present context can 
have a maximal length of Q = k. The rapidities representing possible Q-particle bound 
state are given by 

{ u } = { u + -(Q + l-2j)\j = l,...,Q}, Q = l,...,k. (5.1) 
9 

uw-strings are most naturally classified in as two types; string of length less than k and 
strings of length one with a particular type of complex rapidity. We will refer to these as 
strings with positive and negative parity respectively 18 . The roots constituting possible 
positive parity M\vw string are given by 

{ w } = {v + -(M + l-2j)\j = l,...,M}, 
9 

{ v \ = { v + -(M-2j)\j = l,...,M-l}U{v + -M,v--M\, (5.2) 
9 9 9 

for M = 1, . . . , k — 1 and v real, whereas the negative parity length one string has its center 
on the line ik/g 19 . To each of the v rapidities in the first set of the second line we associate 
two y-roots, one with \y\ > 1 and one with \y\ < 1, while the single y-roots associated to 
v + iM/g and v — iM/g have \y\ < 1 and \y\ > 1 respectively. In terms of the mirror x 
function we have y±M = x(v±iM/g) for positive parity strings, while y\ = x(v+i(k + l)/g) 
and y_i = x ^ v+i (l_ 1 y g ^ ) with v real for negative parity strings. 

18 This is the nomenclature in the XXZ model [17]. 

19 This is not a strange type of solution, merely a consequence of parametrization. In the u-type variables, 
this simply corresponds to negative roots which are perfectly allowed in the Bethe equations. 
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w-strings come in the same types as vw strings, which is no coincidence as we explain 
in the appendix. The set of w-roots making up a M\w string is given by 

{ w } = {w + i(M + 1 - 2j)/g} , j = l,...,M, (5.3) 

where M = 1, . . . , k— 1 and w real or M = 1 with w on the line ik/g for positive respectively 
negative parity strings. 

y-particles are characterized by the property that \y\ = 1. This means their associated 
rapidities v run over the interval (—Ub, In terms of the mirror x function we have 
y~ = x(v) while y + = -t-t, where the ± denotes the sign of the imaginary part of y. 

Negative parity strings turn out to disappear rather naturally from the problem; their 
Y- functions are inverse to those of strings of length k— 1. We will derive this carefully below, 
though this fact already manifests itself at the level of the string hypothesis; by fusing a 
length one negative parity (v)w-stvmg with a length k—1 (v)w string we get a configuration 
with momentum ir, which effectively essentially means the associated Y- functions must be 
inverse in the thermodynamic limit. 

Given a string hypothesis and the corresponding set of Bethe-Yang equations the 
derivation of the TBA equations becomes a text book story. There is a subtlety regarding 
the negative parity states, so for this reason and general completeness we go through 
the details carefully for the interested reader in appendix C. This gives a set of coupled 
integral equations between so-called Y-functions, one for each particle type in the spectrum 
we just discussed; k Yq functions for Q-particles, Y± functions for y particles with positive, 
respectively negative imaginary part, and k Ym( v )w functions for (v)w strings. Of course 
we have a set of Y± and Ym\( v )w functions for both values of a. The index on the Y M u v \ w 
functions runs from zero to k — 1, where we denote negative parity strings separately as the 
zeroth type of strings. However, as indicated above and derived explicitly in the appendix 
(cf. eq. (C.25)) the Y-functions for negative parity strings are inverse to those of length 
k — 1 positive parity ones 

Y 0\(v)w = ^(a) • ( 5 - 4 ) 

k—l\(v)w 

We eliminate YquI\ w from the equations in this fashion, and consequently we introduce the 
following notation. The indices M,N and L run from one to k — 1, while Q,P and R run 
from one to k; repeated indices indicate a standard sum. The * denotes convolution, the 
different types are explicitly defined in eqs (C.10-C.12). This gives the following equations. 

5.2 The canonical TBA equations and free energy 
The TBA equations for Q-particles are 



log Yq = — L8q + log (1 + Yp) * K*% + log ( 1 + ) * (5.5) 



sl(2) B \ 1 y' 

x M\vw 

+ 1'* ( 1 + K%J * + log ( 1 - i ) *Kf 
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where we have an implicit sum over a = 1,2 and (3 = ±. 
The TBA equations for y-particles are 

log ri a) = — log (i + Yq) * k± v + bg - — ^ * + log ^ Vrr- * ■ ( 5 - 6 ) 

The TBA equations for w-strings are 

1^ Y M\ W = log ( : + %) * K-NM + bg (l + rftj * ^-1,M (5.7) 



(«) 



1,A/ 



(5.8) 



+ log i * K M ■ 

The TBA equations for uw-strings are 

^ *ilL = lQ g f 1 + psr) * K-NM + log (l + l£U * IT*. 
1 L_ 

+ log * K M - log (1 + Y Q ) * Kg M . 

1 v («) 

Clearly Yfc_i| ulu and Y k _±\ w enter the above equations in a special way due to eq. (5.4). As 
usual these equations can be simplified and subsequently used to derive the Y-system; this 
requires identities satisfied by the kernels which we show in appendix D.2. 

The free energy of the mirror model allows us to find the energy of the deformed 'string' 
theory as 

E{L) = ~ /^£^log(l + lQ). (5.9) 
J Q=l 

Note that we have only a finite number of TBA equations and a finite number of Y- functions 
contributing to the energy. 



6 Simplified TBA equations and Y-system 

In this section we simplify the canonical TBA equations of the previous section and derive 
the corresponding Y-system. 
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6.1 Simplified TBA equations 
w and vw-strings 

Using (D.26), (D.27) and (D.39) and the canonical equations for log we obtain the 
following simplified TBA equations for uiu-strings 

log = log (1 + Y M+ 

l\vw — l\vw 

) -k S — log (1 + Ym+i) * S 

+ «5 M) ilog ( \ — ^)* s ' 



1-Y+, 

l°g*fc-2|l«l> = log(l + ^-3|t>«;)( 1 + y fc-l|^) 2 *S- ^g (1 +lfc_i)*S, (6.2) 

log Y k _ x \ vw = log (1 + Y k _ 2 \vw) * s ~ log (1 + Y k ) * s , (6.3) 

where Y \ vw = 0. Note how the contribution of Y k _^ vw is still special, but now definitely 
more elegant. For w-strings we similarly find 

f l ~Y-\ 

log Ym[«j = log (1 + Y M+1 \ W )(1 + Y M _ llw ) *s + S M ,i log y - _ j_ J *s , (6.4) 

log>fc-2|«, = ^og(l + Y k _ 3lw )(l + Y k _ llw ) 2 -ks, (6.5) 

log ifc-i| w = log (1 + Y k _ 2 \ w ) -k s , (6.6) 

where Y \ w = 0. 

Q-particles 

Applying (K + 1) _1 defined in (D.29) and using the identities 20 in appendix D.2, we find 
the simplified TBA equation for Y\ 

log Yi = log - j—^- * s - A*s , (6.7) 

1 + Y2 

where 

A = + log (l - ^\ (l - ^\ *K + 2 log(l + Y Q ) * K% (6.8) 
+ log (l - * Km + log (l - y£[ ) * K k -i . 

\ 1 M\vw J V ' 

The kernels appearing in the expression above are defined in appendix D.2. Also here and 
below the sum over a is left implicit. Analogously, for Q = 2, . . . , k — 1, we find that the 
simplified TBA equations are given by 

K 1 + Y Q-l){ 1 + Y Q+l) V y q-i\vwJ 



20 We have assumed identities involving the g-deformed dressing phase to be analogous to the undeformed 
case, but have not proven this. 
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For Q = k we cannot simply apply the usual operator 



(K + l)pl = 6 Pk - s(<5 Pfc _i + 5 Pk+1 ) , (6.10) 

since there is no k + 1st equation. Instead, we apply the operator Sq^ — 25q^-\s. As we 
discuss in appendix E, using special identities and the TBA equations for vw strings we 
obtain 

logY k = 21ogY fe _i*s-log(l + r fc _ 1 )*s + log(\ + ^— *s. (6.11) 

\ k — 1 \vw / 

Naively we might have expected to get something like ( 1 + — ) (l + Y k _i\ vw ) in the 
above equation. Because it is exactly analogous to what happens in the equations for vw 
strings it is very natural that we get a doubled contribution instead, though it arises in the 
far less obvious fashion shown in appendix E. 

6.2 Y-system 
w and vw strings 

Deriving the Y-system for u>-strings from the simplified TBA equation is now a simple 
matter; applying s~ l to the equations above 21 immediately yields 

(\ _y-l\ e ( n "~H) 

= (i + r 2 |«) Wzv^j ( 6 - 12 ) 

Y m\ w Y m\ w = U + W[t»)(l + Y M -i\w) , M = 2, . . . , A; — 3 , (6.13) 
n + -2|^-2 h = (1 + n-8|«)(l + n-i|») 2 , (6-14) 

n-ii^-i k = ( i +^-2iJ- (e-is) 

For -uiw-strings we similarly get 

/1 + Yoi \ /i_y \ e ( u H«l) 



1|ijm> l|t>ui I ^ _|_ J I ^ y 

y jJ|tm; y jfer|wu = ( X + y M+l|t,«;)(l + >Af-l|«^) ( i + ) ' ^' 17 ' ) 

1 fc-2| t ;u- 1 fe-2|^ = C 1 + Y k-3\vw)(l + Yk-l\vwf (j —J , (6-18) 

^fc-llw^/e-llww = ( 1 + 1 fe-2|i ? «;) f j + • (6-19) 

Away from M = 1 we indeed see precisely the XXZ structure advertised in figure 1 in 
the introduction for both the w and vw functions independently. The two Y-functions in 
the lime-green bubble there are Y k _iu v \ w and Y u v \ w and are inverse to each other. We 
implemented their relation at the level of canonical equations already, so that we precisely 

21 As usual / o s _1 (u) = lim[/(u + i/g — ie) + f(u — i/g + it)]. 

e— >0 
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get the doubled contribution of Y^_\u v -\ w on the right hand side of the Y-system for the 
k — 2th Y- function. Note also how the Yq particles are nicely compatible with this structure 
by existing precisely one length beyond the Y M \ VW functions. Dropping the Yq terms in 
the above now clearly gives precisely the structure of figure 2. 

y-particles 

For y-particles only the Y-system for Y_ is of a nice local form. As in the undeformed case 
we can derive them by applying s _1 to the equation for Y_, and subtracting the canonical 
TBA equations for l\w and l\vw strings to get 

1 + Y {a) i 

Y (a) + Y (a)- = Ijtm, __L_ ( } 

- - i+yWl+yi 

l\w 

Q particles 

Taking \u\ < Ub for Q = 1 we have 



Y+Y{ 



n i 



a=l,2 



Y 



(a) 



Y 2 1 + Y 2 

while for Q = 2, . . . , k — 1 the Y-system is 



(6.21) 



Y+YA 



n (i+f^-) 

a=l 2 ^ Q—l\vw / 



Yq+iYq-i (l + Y Q ^)(l + Y Q+l ) 

Finally at Q = k we have 



(6.22) 



y+y- 



n 

2 ^ fc — lj^KJ / 



n 2 -i 1 + Yk- 



(6.23) 



Taking the above equations together, we have the full the Y-system for the quantum de- 
formed mirror model illustrated in figure 4. The four Y-functions coupling to Y& at the top 
of the diagram are precisely the four Y k _^ vw terms - two for each a - in the above equation. 
As in the undeformed case the Y-system admits analytic continuation onto the complex 
u-plane with cuts where the TBA equations unambiguously fix the corresponding jump 
discontinuities of the Y-functions. Both the TBA equations and the Y-system presented 
here manifestly tend to their undeformed cousins in the limit k — > oo. 

7 Ground state energy 

In this section we will find the constant solution of the ground state TBA equations and 
show that the ground state has zero energy. To do so we follow [68] where this was shown 
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for the undeformed ground state. The idea is that the ground state equations are solved 
by 

Y Q = 0, Y ± = l, (7.1) 

with constant Ym\( v )w functions. However, the TBA equations for Q particles are singular 
at the point Yq = and so we need to regularize the solution and take a limit to obtain 
our result. To this end we write 

Yi a) = 1 + hA^ + .... (7.2) 

As in the undeformed case, the small h expansion of the other equations is consistent with 
Y+ = Y- and constant Ym\( v ) w functions. Then the ground state TBA equation for Q 
particles, eq. (5.5), shows that 

log Y Q = -2 log h * (K yQ + Kl Q ) + finite terms . (7.3) 

Since K_ + is normalized to one, this shows that Yq ~ h 2 BQ where Bq is left 
to be determined. We then need to solve the TBA equations for vw and the completely 
decoupled w strings with constant Y-functions, and use the resulting vw functions to fix 
Bq. The equations for w and vw strings take the same form, namely for Y M u v \ w = Am 
the simplified TBA equations reduce to the Y-system 

A\ { = (1 + A M + A M +1 ) , (7.4) 
A\_ 2 = (1 + A fc _ 3 )(l + A k ^) 2 , (7.5) 
A\_ 1 = {l + A k . 2 ), (7.6) 

with Aq = 0, which in principle has many solutions. However, requiring the solution of 
this Y-system to solve the canonical TBA equations singles out the solution 

Am = M(M + 2) , (7.7) 
A k _ x = k - 1 , (7.8) 

as can be verified by explicit computation. To determine Bq we need to solve 

log B Q = -L£ Q + 2 log (l + ^) * + 2 log (1 + A k ^) * K k v ~l' Q ■ (7.9) 
Noting that is normalized to one, we easily find 

k-2 

S.2 



M=l 



log B k = -L£ k + 2 log f 1 + m(m+2) ) + 2 lo S lh = ~ L£ k + 2 log 2k . (7.10) 



Rather than computing the other terms explicitly, we can use the small h expansion of the 
Y-system (6.21-6.23) to conclude 

Y Q ~ 4h 2 Q 2 e- L£ Q , for Q = l,...,k. (7.11) 
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The associated energy is then given by 

k 



=1 



E~Ah 2 Y,Q 2 I ^e- L£ ^ . (7.12) 



2tt 



Since we sum a finite number of convergent integrals, in the limit h — > the ground state 
energy is manifestly zero 22 . 



8 Conclusion 

In the present paper we have derived the ground state TBA equations of the quantum 
deformed AdS^ x S 5 mirror model, based on a psu ? (2|2) 2 invariant S-matrix. As we have 
seen, these equations show a rich structure, summarized in figure 4. The next natural step 
is to generalize these equations to excited states along the lines of [25-29] , and as indicated 
we will provide the necessary tools to do so - an asymptotic solution - in an upcoming 
publication [55]. 

The main important feature of the equations we have found is that there is a finite 
number of them. This will hopefully allow us to gain further valuable insight into the 
AdSs x S 5 mirror model by considering our equations in an appropriate limit. Having 
only a finite number of TBA equations may help avoid technical obstacles which arise 
in the strict undeformed limit. As one potential example, when the above equations are 
generalized to excited states an obvious question will be that of wrapping corrections at 
small coupling, for example for the deformed analogue of the Konishi operator. With a 
finite number of TBA equations, it might be possible to obtain an analytic expression for 
wrapping corrections at a fixed value of k 23 without relying on the clever but unfortunately 
limited approach of [30, 35] which appears to work specifically just beyond first wrapping in 
the s[(2) sector. The idea would then be to find a sequence of wrapping corrections at fixed 
values of k, and hope that this naturally extrapolates to the undeformed case, giving higher 
order wrapping corrections for SYM operators in the limit k — > oo, previously unaccessible 
analytically. 

Coming to more general questions, as we saw above the quantum deformed model at 
roots of unity shows qualitatively rather dramatic differences to the undeformed case. In 
contrast we can expect deformations with real q to be 'smoother' and can ask whether 
the resulting model has interesting properties. Considering the amount of material already 
presented, we decided against trying to incorporate this in the present paper. However, 
we should note that the answer to the above question is positive, though admittedly less 
dramatic than was the case here. We would like to briefly elaborate on this. The XXZ 
model considered as a real deformation of the XXX model has | A| > 1, and in this case the 
TBA equations [18, 62] are simply modified by being appropriately defined on an interval 

22 The divergence of this expression for L — 2 observed in the undeformed case [68] does not occur here 
since we only have a finite number of integrals. The sum does grow with k however, and for L = 2 the limit 
k — > oo and h — ¥ do not commute. 

23 The previous obstacle of inversion of an infinite dimensional matrix of integral operators now at least 
becomes finite dimensional. 
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rather than a whole line, while otherwise their structure remains intact. Of course a similar 
thing happens in our model when q is real, but the fact that we have an additional coupling 
constant in the game results in something special. Firstly, for fixed q there is a natural 
bound to the coupling constant g above which the model is no longer well defined, and 
secondly the support of Y± functions varies with respect to that of the other Y-functions 
and becomes equal to it exactly on the boundary for g, which seems to result in a Y-system 
without discontinuities provided this limit is well defined. These comments of course also 
apply to the deformed Hubbard model subsystem. We hope to address this in more detail 
in the near future. 

Also, let us note that the thermodynamics of the deformed Hubbard model itself might 
have interesting features from a condensed matter point of view, and apparently provides 
the next simplest setting beyond the XXZ model to study the TBA of quantum deformed 
systems. It would be interesting to look into the thermodynamics of the quantum deformed 
Hubbard and associated models 24 in more detail along the lines of e.g. [69]. 
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A The S-matrix 

In this appendix we summarize the properties of the g-deformed S-matrix. We use E^ to 
denote the 4x4 matrix unity, i.e. a matrix with a one in the (i,j)ih entry and zeroes 
everywhere else. Next, we introduce the following definition 

E kUj = {-l)< l >^E ki ® En , (A.l) 

where e(i) denotes the parity of the index, equal to for i = 1,2 (bosons) and to 1 for 
i = 3, 4 (fermions). The matrices E^uj can be used to write down invariance with respect 

24 While in undeformed Hubbard-like models we can directly go from the 'mirror' type model (Hub- 
bard) to 'string' type models (corresponding to noncompact support for the Y± functions in the standard 
parametrization) it appears that with the root of unity deformation the string hypothesis for the 'stringy' 
type deformed Hubbard model requires more than the obvious modifications. 
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to the action of two copies of U q (su(2)). If we introduce 

Al = #1111 + |#H22 + ^(2 ~~ <? 2 )#1221 + 2^2112 + |-E^2211 + #2222 , 

A2 = --E1122 — 1^1221 ~~ 2g^ 2112 ~*~ 2^ 2211 ' 

A3 = #3333 + |#3344 + -(2 - q 2 )Es443 + -#4334 + |#4433 + #4444 , 

A4 = -^3344 — ^#3443 — 2^^ 4334 2^ 4433 ' 

A5 = £ll33 + #1144 + #2233 + #2244 , (A-2) 

A6 = £3311 + #3322 + #4411 + #4422 , 

A7 = -©1324 — <7#1423 #2314 + #2413 , 

Q 

A8 = #3142 — <?#3214 #4132 + #4231 , 

q 

Ag = #1331 + #1441 + #2332 + #2442 , 
A10 = #3113 + #3223 + #4114 + #4224 , 

the S-matrix of the g-deformed model is given by 

10 

Si2{pi,P2) = ^2a k (p 1 ,p 2 )A k , (A.3) 
k=i 

where the coefficients are 
01 = 1 , 

_ _^ _j_ 2 x 1 (1 — x 2 x~\)(x~{ — x 2 ) 



as 



qxf(l-x 1 x 2 )(x 1 -x%) 
U2V2 x+ - x 2 
U1V1 x~[ — x 2 



U 2 V 2 xt-x 2 2U 2 V 2 x 2 (x+ -x+)(l 
a 4 = -q——^z + 



UiVix 1 — x 2 qUiVix 2 (x 1 — x 2 )(1 — x 1 x 2 ) 

05 = X '7 X * TT (A.4) 

^UiVi{x 1 -xj) 



a 6 



a 7 



ag 



aio 



Vqu 2 v 2 (x^ 



X q 



x 1 x 2 

ig (xf - x^)(xf - x%)(x% - x 2 ) 
2 y/qUiVi{x^ -x^")7i72 
2j [/2V2 xj" 3 ^ - xj")7i72 
5 qixfxj (xy - x 2 h )(x^x 2 - 1) 
(xj~ - 3^)72 
{xl - x+ ) 7 i 

^"2^2(X2 ~X 2 f )7l 
flVl (Xj - X^)72 
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Here the central charges are given by 



u 



Vi 



, x t x i + £ 
'—+11' 



(A.5) 



and the parameters 7, are 



(A.6) 



The dependence of the S-matrix on the variables 7$, i = 1,2, is gauge-like. Indeed, intro- 
ducing the diagonal matrix Tj = diag(l, 1, 7i, 7»), we find 



rr 1 ® r 2 - 



^12(^1,^2) 



(A.7) 



where SJ2 1 is the S-matrix where 71 and 72 are set to one. 



Let us summarize the most important properties of the S-matrix. The S-matrix satisfies 

• The Yang-Baxter equation; 

• The unitarity condition 

#21(22, z l )Si 2 (z 1 , z 2 ) = 1; (A.8) 



• The transposition property 

s t {z u z 2 ) = t g n S(z lt z 2 ) n-Ha , 

where I 9 = {—l) ei ^E\ ® E 3 - is the graded identity and Q is given by 
n = exp — \Ei ® + F 1 ® + E 3 ® F 3 + F 3 ® E 3 



(A.9) 



(A.10) 



Here Ei, Fi are positive and negative roots of two su(2)'s. The graded identity com- 
mutes with S7. In the limit q — > 1 one finds 



lim 



Q,S(zi,z 2 )n = lim S(z\, z 2 ) . 



Note also the formula x + (z, q) = —x (—z, l/q). 
Generalized physical pseudo-unitarity 



where 25 



S( Zl ,z 2 ? = BS- 1 (zl,z*)B- 1 
B = A® A, with A = diag(ai, 01) 



(A.ll) 
(A.12) 



3 Our matrix B is Hermitian as is required for pseudo-unitarity [56]. 
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and o\ 




is the first Pauli matrix. 



In addition, the S-matrix has a unitary spectrum on the real line of the string and 
mirror theory, but is generically not unitarizable by a local basis transformation. 
In fact the many body S-matrix on the real string line has non-unitary eigenvalues. 
Nonetheless, on the mirror line the many body S-matrix appears to have unitary 
eigenvalues always 26 . This is equivalent to being quasi-unitary, where quasi- unitarity 
means that the matrix B above is positive definite; B = OO^ [57]. However, there 
does not appear to be a matrix B of this form which is also factorizable over the 
one-particle basis, and as such we have not rigorously proved quasi- unitarity on the 
mirror line. 

For coincident arguments the S-matrix reduces to the permutation. 

It is compatible with crossing symmetry. If we introduce the following charge conju- 
gation matrix 



C = 



I 





— iq 


iq 


-1/2 













\ 









1/2 






q 
-q- 1 ' 2 / 






,1/2 



the crossing relation reads 27 



Su(zi, z 2 )C 2 S t 1 2 2 (z 1 ,z 2 - w 2 )C 2 1 



?(xr-xj)(l 



(A.13) 



The complete S-matrix comprising two copies of S can be written in the form 

, _ ,i l 

s = s M2) s®s, s sm = ^-—^-l^ ■ x \ _ x 2 - ! _ _V > ( A - 14 ) 

where g) stands for the graded tensor product and a is the dressing phase [63, 64]. Substi- 
tuting this representation for S into the crossing equation [65] , we deduce that the dressing 
phase must obey the following equation 

a(z 1 ,z 2 )a(z 1 ,z 2 -u 2 )=q- 1 ^ Xl _ ~ 4 f . (A.15) 

26 We have checked this numerically for the first few n-body S-matrices. 

27 Under the crossing transformation x — > l/x so that the central charges transform as U 2 — ¥ l/U 2 
and V 2 ^l/V 2 . 

In [12 the factor L + in the S-matrix was replaced by =jt = -4: i — . This is a minor change which 

leads to the corresponding modification of the crossing equation for a. 
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Using the unitarity relation the previous formula also implies that 



a(z 1 + u 2 ,z 2 )a(z 1 ,z 2 ) = g' 1 ^ • ^ ^- . (A.16) 

x 2 x 1 x 2 1 + 

The formulae (A. 15) and (A.16) are ^-deformed analogue of eqs. (2.8), (2.9) and (2.10) 
from [66]. A natural solution of this equation has been obtained in [12], where an explicit 
expression for a can be found, generalizing the DHM representation of the undeformed 
dressing phase [67]. 

The dressing factor a p Q which describes scattering of P- and Q-particle bound states 
can be obtained from the dressing factor for fundamental particles by means of fusion. The 
corresponding crossing equations have almost the same form as in the undeformed case, cf. 
eq.(2.14) in [66], 



h PQ ( Zl ,z 2 ), (A.17) 



a p Q(z 1 ,z 2 )a p Q(z 1 + u, 2 , z 2 ) = <T P<? P h p ^(z 1 ,z 2 ) . (A.18) 



Here x\ and x 2 solve the bound state conditions 



?(^ + ^)-' p (^ + ^)=(' F -?) H 
^^(^H^H 1 ' (A - 20) 

and we have introduced the crossing function 

i 1 p-i 

2_ 



foPQ _ ^2 



II S Q-p+V . (A.21) 

+ — -_1 



where <Sq is defined in (D.2). 



B String hypothesis for the deformed Hubbard model 

As explained in the main text, in the limit R — > oo we obtain Q-particle bound states. 
As we will show here, going down one level and considering the limit K l — > oo yields vw- 
strings, while going down two levels and considering the limit K 11 — > oo gives w-strings. 
While not in the strictly logical order, we will start by discussing the string hypothesis for 
w-strings because it gives the clearest illustration of the effects of the deformation. 

B.l y-particles, u>-strings, and the XXZ model 

To describe the string hypothesis at the level of the second auxiliary Bethe equation, we 
need to know which type of y-roots can appear there. This is derived in the section below, 
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where we will see that \y\ ^ 1 necessarily leads to the emergence of a vw-sthng so that y- 
particles are characterized by the condition \y\ = 1. With this condition, the equations for 
w particles become nothing but the Bethe equations of the inhomogeneous XXZ spin chain, 
where the inhomogeneities are unitary defects. As such the results here closely mimic the 
string hypothesis for the XXZ spin chain [17, 18], but the derivation is slightly different in 
form. 

The discussion is simplest if we consider the second auxiliary Bethe equation in the 
form (3.2) 



K u K U1 2 



w k - qv o . 
i=i i=i 



(B.l) 



where in this notation \yi\ = 1 implies that the Vj are positive. We want to consider the 
types of solutions these equations admit in the limit K 11 — > oo. Now for q on the upper 
half unit circle and v positive we have 



qw — v 



w — qv 



> 1 for Imw > 



(B.2) 



and vice versa. Thus, the product over y-particles diverges or goes to zero unless w is real. 
Without loss of generality we will start with a w root, say wi, for which the product in the 
bracket below diverges 



n 



wi - qvt 



K m 2 
J-j- VTj - </-\Y, 



i=l 



'Wi — Wi 



(B.3) 



The only way this root can be part of a solution of the Bethe equations is if this divergence 
is compensated by a zero in the second product in the Bethe equation, say 

wi - g 2 w 2 = . (B.4) 

However we should now also ensure that there are no problems in the equation for w 2 



1 



n 

,1=1 



qw 2 - Vi 
w 2 - qvi 



K m 2 
T-r w 2 - q^Wi 

11 q 2 W2 _ w . 



(B.5) 



This can be determined by multiplying the two equations (B.3) and (B.5), giving 

K u K IU 2 2 

t-t qwi - Vj qw 2 - Vj j-i- wi — q Wj w 2 - g z w« 
Al wi - qvi w 2 - qvi q 2 Wl _ w . q 2 W2 _ w . ' 



(B.6) 



whereby the divergent term in (B.5) originating from (B.4) cancels out. At this point there 
are two possibilities. Either the equation for w 2 is finite in the limit K 11 — > oo, which is 
equivalent to 

I , (B.7) 



qwi 


— v qw2 — v 




q 3 W2 — v 


Wi - 


qv w 2 — qv 




w 2 — qv 
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giving a string of length two, or this product is greater than one and the process continues 
to give a longer string. In the latter case we necessarily have to involve another root to 
find 

w 2 -g 2 w 3 = 0, (B.8) 

and analyze its equation in turn. Continuing along these lines we see that a general string 
configurations of w-roots are given by 

{w} = { 9 M+1 - 2 M, j = l,-..,M. (B.9) 

where M denotes the length of the string. Consistency of the resulting string solution 
requires that the center of the string is real, as can be seen via fusion, see appendix D.l. 
This is not the only constraint however. 

As we saw above, in order for our string configuration to keep growing, we must keep 
a divergent product at every stage of the derivation. This means that in order to obtain a 
string of length M we must have 



n 

k=l 





— V 




qi(v-q M w) 


Wfc - 


qv 




q2j v _ qM w 



> 1, for j = 1,...,M - 1. (B.IO) 



Taking into account that v is positive and q = e l fc , we find that the last condition is 
equivalent to the following requirement 

wcos(M^) < wcos ((AT - 2j)^) , for j = 1, M — 1 . (B.ll) 

Given these conditions, we can have strings of length one through k — 1 for a positive 
center (w > 0) and a single type of string of length one with negative center (w < 0). 
We will refer to the positivity or negativity of the string center as positive and negative 
parity respectively. In terms of the u-plane parametrization, positive parity strings have 
real rapidity w, while the rapidities of negative parity strings lie on the line +ik/g. 

Note that the positive parity string of length k lies just on the edge of these conditions 
and is not allowed; there would have been no divergence to begin with. However, were we 
to perturb slightly away from k being an integer and let k be m + e, of course it could 
become an allowed string configuration; correspondingly in the limit where k is exactly an 
integer this string configuration has fixed momentum tt and trivial scattering with all other 
roots, which is another way to see it plays no role in thermodynamics. 

The constraints on the length of strings in the XXZ model were first postulated in 
[17] and later shown to be equivalent to normalizability of the Bethe wave function of the 
associated string configuration [70] . The argument of normalizability of the wave function 
would of course also apply here. 

B.2 uw-strings 

Now that we have seen how to analyze the string solutions of the equations for w roots, let 
us take a look at the equations for y roots in the limit K l — > oo. To analyze this limit, we 
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need to analyze the product over Q particle excitations in (3.1). This product consists of 
terms of the form 

S Qy = q Q/2 x LzA x f x L, (b.12) 
x + — y V x 

where x describes a Q-particle bound state with real rapidity u. Due to the conjugation 
properties in the mirror theory (4.12) the ratio x + /x~ is real. Hence the modulus of S® v 
is given by 

\ s ^ = y*-^y-\ x -Jy\ 2 £\ l -y* x+ -^. (B .i3) 

1 1 y*-^y- X + X - \ y 2\ + 2±_ y * x+ _^ 

Hence this modulus is smaller than one provided 

1-— Vl-|y| 2 )<0. (B.14) 
x J 

Now for any bound state number Q = 1, ... ,2k — 1 and real u we have x + /x~ < 1. So 
\S Q y\ < 1 for \y\ > 1 and vice versa. As in the undeformed case, this implies that in the 
limit K l —7- oo the product over the momentum carrying particles in (3.1) diverges when 
\yk\ < 1 and goes to zero when \y k \ > 1. Auxiliary roots that have \yk\ = 1 are identified 
as y-particles and w-strings exist in their presence as we discussed above; here we assume 
that \y k \ ^ 1. 

Without loss of generality we take \yi\ < 1. Then in the limit K l -4 oo eq.(3.1) 
becomes 

1 *L sinh |f (v k -Wi- I) 
l = -xTT == 9 —, B.15) 

which requires one of the numerators in the product on the right hand side of the equation 
to vanish. This can be achieved by taking for instance 

«i-tOi-- = 0. (B.16) 
9 

We should also satisfy the equation for w\ however, which now gives 

1 = x n ¥k(Wl ~ Vl + ^ ff ¥k(Wl ~ ^ ~ 1 } , (B-17) 
i =2 sinh K^i - v i ~ |) ~=i sinh (wi -Wj + f)' 

so that we should take 29 

lBl-V2-- = 0. (B.18) 
9 

Now 2/2 has become involved in the game, and there are two options. If 1 2/2 1 > 1 the 
equation for yi is 'satisfied' and we get a l\vw string. 



Taking a denominator in the product over w roots to vanish would give vi + - = W2 which contradicts 
the original considerations in (B.15), since the zero in the equation for wi by assumptions scales as the one 
in (B.15). 
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In case \y 2 \ < 1 we get further conditions which we explicitly work out one further 
level before giving the general pattern. Assuming 1 2/2 1 < 1> analogous to (B.15) and (B.16) 
we take 

v 2 -w 2 -- = 0. (B.19) 
9 

However we now have the equation for w 2 as well, which tells us 

1 = x TT Sinh ™ {W2 ~ Vi + 9 ) TT sinh 1^ - Wj - |) ^ (B 2Q) 

~T 2 sinh If ( W2 ~ Vi ~ 9 ) /=i sinh - ^ + f ) ' 

so that u>2 — — ^ = 0. However we should not forget that from (B.18) and (B.19) we 
now also get 

2i 

w 2 = w 1 , B.21 

9 

meaning we should additionally have a numerator vanish 30 by imposing a condition 

w 2 - v 4 + - = . (B.22) 
9 

At this point we have the configuration 

2i % i 1% 

v\ = v 2 = w\ = w 2 + - = u 3 H = V4 , (B.23) 

9 9 9 9 

The configuration (B.23) forms a 2\vw string if 1 2/3,4 1 < 1. If 1 2/3 1 > 1 or \y^\ > 1 we get 
more conditions, and in general we can end up with an M\vw string given by the following 
configuration of roots 

{ w \ = iv + (M-l)-,...,v-(M-l)- 
9 9 



{v} 



v + M-, v + (M -2)-,...,v- (Af-2)-, v-M-) . (B.24) 
9 9 9 9 



M-l 

Note that there are two y- roots yj,l/ Vj associated to each of the underbraced rapidities 

Vj = v + (M-2j), j = l,...,M-l 

in a M\vw-strmg. We denote by yu and y~M the roots associated with the first and 
the last rapidities in (B.24), respectively. According to our construction, \yu\ < 1 and 
\y~M\ > 1. In total there are 2M y-roots. 

As for w-strings, the reality conditions for the string centers follow by fusing the S- 
matrices and insisting that the S-matrix for the total M\vw string is unitary. Fusing the 



30 

denominator 



Note that we have W2 = Wi — ^ only 'after' imposing (B.19); it cannot fix the original vanishing 
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S-matrix S® y over the constituents ?/-roots of an M\vw string gives 

TT sQy(u, v) = qQ M \ ~ VM \ ~ y ~ M x - TT *--Vi x ~-k *l 

ieuU " X ~ VM X ~ y - M X ~ ti X ~ ViX+ ~ i X 

M-l - , J_ _ _ J_ 
QM x ~VMX - y-M x TT x ^ x- j» V i 

<r 11 . ! „, ' 

Vi 



X 


- UM X 


- y- 


_i_ 

-M % 


X 


yM x 


7/ 

y- 


-M x 


x~ 


- VM X~ 


- y- 


-M X + 


x+ 


~ VM 


- y- 


-M X~ 


x~ 


~ VM X~ 


- y- 


-M X + 


x+ 


~ VM X+ 


-y- 


■M X~ 


x~ 


~ VM X~ 


- y- 


-M X + 



=T x + x+ ~ Vi 



qm x -yux -y- M x' yr q k - q m a e k 

" q^e k — q M ~ 2l e * 
M-l 

QM x — VM x — y-M x -p-i- M-Q-2i a I \ m 

(T — r \ \q v S Q+ M-2i{u-v) B.25 



y M x - y-M x n 
VM x+ - y-M x~ 



i=i 
M-l 



11 S Q+M - 2i {u-v). (B.26) 



Requiring this S-matrix to be unitary requires that y* M = — — and n to be real modulo 
ik/g. 

At this point we do not have a restriction on the string length yet. This is because at 
every stage there are two y roots associated to each rapidity, which means we can always 
choose the one with the desired properties. However, going beyond length k it is not hard 
to see that we must necessarily have coincident roots. For w roots however, we had no 
length k string either. Again here the peculiarities of y-roots do not allow us to dismiss 
the length k string; we can obtain a bona-fide divergence and repeat the arguments above. 
However, similarly to the case of w-strings the would-be k\vw string has fixed momentum 
and scatters trivially with everything else. This means our strings can have length up to 
k-1. 

Next, coming back to the location of the center, it is easy to see that by translating a 
given root configuration by 2ik/g the set of all but the outer two y- roots is immediately 
invariant. Now since we have the constraints \vm\ < 1 an d y* M = a ^ so ^ e outermost 
y roots cannot change either, and hence the center of the string can be taken to lie either 
on the real line or on the line ik/g. 

At this point it appears we can have negative parity vw strings of length up to k — 1, 
while for w strings this was not the case. Also this strange fact arises due to the curious 
property of the y-roots mentioned earlier; we can naively choose the correct y roots to 
make our string hypothesis work. However, the root configurations that lead to negative 
parity strings of length greater than one are not consistent at the at the level of w roots, 
where they violate the conditions we derived above 31 . This in particular implies that the 
Bethe wavefunction associated to a negative parity vw string of length greater than one in 
the deformed Hubbard model subsystem would not be normalizable. Hence we are led to 
exclude negative parity string of length greater than one, and end up with a string content 
exactly analogous to the second auxiliary level. 



31 Of course the conditions on the w roots immediately imply the conditions on the vw strings derived 
above independently. 
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Parametrization of y-roots 

As mentioned earlier, x{u) runs over the lower half unit circle as u runs from — Ub to u^. 
Therefore, a y-particle with negative imaginary part is parametrized by x(u), while a y- 
particle with positive imaginary part is parametrized by ^hjx, with — uj, < « < u&. In other 
words 

V + (v) = ^, y'(v) = x(v). (B.27) 

X(V) 

Next, in an M\vw string of positive parity the outermost y roots are y±M, associated to 
the rapidity v±iM/g, with the requirement that \hm\ < 1 an d \V-m\ > 1- These roots are 
therefore parametrized as 

y+ M = x(v ± iM/g) . (B.28) 

From the properties of x(u), see e.g. figure 7, it follows that Imy»v < for any M = 
1, . . . , k — 1. The negative parity string of length one has two y-roots associated to v + 
i(k± l)/g, with < 1 and \y~i\ > 1. These roots are therefore parametrized as 

* = X(U + + ' yZl = x(H»(l-l)/ g ) ' (R29) 
Here we have Imy^ > 0. 



C Derivation of the TBA equations 

The Bethe-Yang equations in the thermodynamic limit 

As indicated in the main text, given a string hypothesis and the corresponding set of 
Bethe-Yang equations the derivation of the TBA equations becomes a text book story. 
For completeness we present it in this appendix. The first step is to find the integral 
equations for the hole and particle densities, the analogue of the Bethe-Yang equations in 
the thermodynamic limit. In the present context there is a subtlety regarding the definition 
of counting functions, so for this reason as well as general continuity let us quickly go 
through the general story. 

The string hypothesis classifies the types of solutions that can arise to the Bethe-Yang 
equations in the thermodynamic limit. This means that the products over elementary 
excitations in the Bethe-Yang equations should arrange themselves as products running 
over the allowed string complexes containing the relevant excitation, where each string 
type of complex can occur a given number of times. In other words we have a product over 
the (r) allowed types of strings labeled by M, the number of strings of type denoted Nm, 
and a product over the constituents of individual strings, that is 

K r N M 

ij «<)-»• n n n ^>o- ^ 

i=l M=lm=ljeM m 

Due to the structure of the string complexes, the S-matrices involved in the products over 
individual strings fuse to give the S-matrix that describes scattering of the entire complex 
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m one go 

H S{u uUj ) -> S M \u u u m ) , (C.2) 

where we labeled the central rapidity of the string M m by u m . This gives us the Bethe-Yang 
equations for scattering of elementary excitations with string configurations. To describe 
scattering of string complexes between each other, we take a product of these Bethe-Yang 
equations over the external elementary excitations making up a string of length N 

r N M r N M 

e ~ ipiL = n n sM ( u h u m) -> z~ imL = n n s NM ( UN , Um ), (c.3) 

M=lra=l M=lrn=l 

where we labeled the momentum and rapidity of the string as pjy and ujy respectively. For 
the auxiliary equations, the analogue of e m is the S-matrix for scattering of the higher 
level root with the root under consideration. The fused S'-matrices implicitly used here are 
explicitly given in the next appendix, section D.l. 

Next, we take the logarithm of the Bethe-Yang equations. Through a fixed choice of 
branch for the logarithm this introduces an integer / in each equation which labels the 
possible solutions of the equation 

r N M 

2nI N = L PN -iY,Yl lo § SNM ( U N, um) , (C.4) 

M=\m=\ 

Now in the thermodynamic limit the solutions of the Bethe-Yang equations become dense 
and it makes sense to generalize the integer 7 to a function of the relevant momentum, the 
counting function, interpolating between the integers corresponding to different solutions 
of the Bethe-equations. It is important that this counting function is monotonic so that if 
increasing, the derivative of this counting function gives us the density of possible solutions 
as a function of momentum. If the derivative is negative, since we would like to think of 
densities as positive, we should change the sign of the definition of the counting function. 
By definition, the total solution density is then a sum of the particle and hole densities, 
corresponding to solutions occupied and left vacant respectively. 

Simultaneously, on the right hand side of the logarithmic Bethe-Yang equations the 
sums over excitations solving the Bethe equations turn into integrals over particle densities 

r N M r 

E E -+ E / ( c - 5 ) 

M=l m=l M=l J 

In short, we take the logarithmic derivative of the Bethe equations and equate the result 
to the sum of the particle and hole density for the excitations under consideration. This 
gives us integral equations for particle and hole densities; the thermodynamic analogue of 
the Bethe-Yang equations 

sign (^) ( PN ( v ) + p N (v)) = ^ d ^+f] [ duK NM (v,u)p M (u) , (C.6) 



7 M=l 
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where K NM (v,u) = log S NM (v,u). The above sign is precisely such that we have a 

monotonically increasing counting function in accordance with the string hypothesis, see 
e.g. [60]. This sign is especially important in the present context; it comes in due to 
remarkably nice properties of the scattering kernels following from the special nature of 
the deformation. Namely for w and vw strings we have 

^lo g S^(u,v) = -^lo g S^(u,v), (C.7) 

^logSfW) = ~logS^-\u,v), (C.8) 

^logS 00 (u,v) = A log 5 11 («,«), (C.9) 
au au 

where the last equation holds simply by the definition of type w and vw strings, and 
X is a generic label denoting any relevant S-matrix. These properties immediately show 
that negative parity strings have counting functions defined with a relative minus sign with 
respect to their positive parity cousins. 

In light of the above properties it should be clear that the canonically defined kernels 
for negative parity w and vw strings are given by minus the kernels for their length k — 1 
positive parity cousins. Hence we introduce the following temporary notation: repeated 
indices of type M (M, N and possibly L) do not quite denote the standard summation, 
but rather a M bu = Ylp=i aP ^p ~ a °^k-i- Repeated indices of type Q (Q, P and possibly 
R) denote a standard sum with Q running from one to k. Let us also define three types of 
convolution 

/oo 
dtf(u,t)h(t,v), (CIO) 
-oo 

dtf(u,t)h(t,v), (C.ll) 

/-U b j-OO 
dtf(u,t)h(t,v)+ / dtf(u,t)h(t,v). (C.12) 
-oo J u b 

We mention again that all S'-matrices are derived and defined in section D.l. With these 
definitions we explicitly have the following equations 

Q-particles 

PQ + ^ = ^f u +K %)*P? + t (E^^ + ^^SU) ■ (C13) 

a=l \/3=± / 

where generically we have K A (u,v) = log S A (u, v), except of course for K® V M at 

M = 0. 

y-particles 

T (P^ + = * PQ + K AI * P W W + K N * P W , (C.14) 

where Kf is special; Kf(u,v) = ^ t £logS^ y (v,u). 
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w-strings 



P% + til = + P { - ] ) + K*» * P W 

(a) _(a) _ (a) . -(a) 
/'olio + Po\w ~ Pk-l\w + Pk-l\iv 



M = h 



,k-l, 



(C.15) 
(C.16) 



where we made use of (C.7) to rewrite the kernels in the second equation, and of the fact 
that the counting function for a Q\w string is defined oppositely from those of the other 
k — 1 to ensure positive particle densities. 



vu>-strings 

(a) (a) _ rM Q 



K M i{pf+p^)-K MN -- {a) 



(a) -(a) _ (a) -(a) 
Pq\vw Pq\vw Pk— l\vw Pk—l\vw 



Pn\ vw > M = l,...,k-1, 

(C.17) 



log Sxv M {v, u). Also here the counting functions for strings of 



where K^l^(u,v) 
type are opposite to those of the others. 

Free energy and TBA equations 

To derive the TBA equations using the above integral equations for densities, we need to 
minimize the free energy at temperature T = 1/L given by 



du 



2fc-l 



£ Q iu)p Q {u) - J £(-1)^(11) + p#(u)) 



(<*)/ 



a=l 



S 

L 



First we write the free energy compactly as 



T y (L) = / du £ 



£iP. 



L 



(C.18) 



(C.19) 



where £j and ■yj are nonzero only for Q- and y-particles respectively. Next, varying the 
Bethe equations (C.13)-(C.17) gives 



5pi(u) + 5pi(u) = Kij -k Spj , 



(C.20) 



where we note that the implicit sum over j can be of the special kind introduced above. 
Let us explicitly mention that for w and vw strings we have the special case 

fiPo\(v)w( u ) + s Po\(v)w( u ) = S Pk-l\(v)w( u ) + s Pk-l\(v)w( u ) = Kk-lj*$Pj , (C.21) 

which we will come back to shortly. Varying the entropy function then gives 

5s( Pl ) = (ei - ni)S Pi + log (1 + e^- e ') K ti * 6 Pj , (C.22) 

where we note again that the sum over j can be special, while in the sum over i we need to 
take into account (C.21) but that aside from this it is a regular sum. The pseudo-energies 
e,- are defined as 



Pi 
Pi 



(C.23) 
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Upon imposing 5F^{L) = 0, we get the TBA equations 

ej = L£j- log (1 + e i7 *- £l ) * Kij . (C.24) 

For Q\w or Q\vw strings, we should take into account the special sum we introduced above. 
This means the concise notation of Kiq actually stands for —Ki k—i, and since w and vw 
strings carry no energy 32 , we then immediately see that 

e 0\(v)w = — e k-l\{v)w (C.25) 

so that we have just k — 1 independent pseudo-energies to consider for w and vw strings 33 . 
Note also that e y ± is defined only for it G (— Ub,Ub). 

Taking into account that the entropy function can be written in the form 
s(fH) = ~~ log (1 + e^) + (ei - i^pi + log (1 + e**-*) K %3 * Pj , (C.26) 
at the extremum (C.24) the free energy is given by 

^ (i ) = -*/d«£iflo g (l + e<-.). (C27) 

i 

Finally, since the ground state energy of the g-deformed 'string' theory is related to the 
free energy of the q-deformed mirror model as 

E 1 {L)= lim ^(L), (C.28) 

R—>-oo K 

we get 

E 7 (L) = - [ duj2~^ g (l + e-^) . (C.29) 
J Q=l 

The explicit form of the TBA equations (C.24) depends on the particle type; we give 
them in full detail in the main text. There we put them in a form where we have introduced 
^-functions as 

= Y Q , (C.30) 
e «MlW« = Y M]{V)W , (C.31) 
= Y± , (C.32) 

and note that 7 = m for y-particles [4]. We would like to emphasize that the k — 1 w and 
vw Y-functions enter the equations in a special way because of (C.25). From this point 
onward in the main text, the indices M,N and L run from one to k — 1, while Q,P and R 
run from one to k, and repeated indices always indicate a standard sum. 



32 0r viewed differently, their 'energies' are related via (C.7-C.9). 

33 Their precise relationship is different if we were to include chemical potentials, we may come back to 
this in more detail in our upcoming publication [55]. 
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D S-matrices and kernels 



D.l Fusion of S-matrices 

In what follows, unless otherwise indicated the indices M, N and L run from one to k — 1 
while the indices Q,P and R run from one to k. 

S M and S M n 

We will encounter the fusion of the basic S-matrix in the equations for y-particles, w, and 
vw strings, so let us discuss it first. The basic S-matrix Si 



sinh tS~{u — v — i/q) 

Sl (u-v) = —%) -4f{, D.l 

sum \u — v + i/g) 

fuses as 



Sm(u-v)= {I 5x^-1,)= 11 Si(u-v 3 ) = 

ie(v)w je(v)w 2kK /y/ 



(D.2) 



where the product is over the w roots making up an M\(v)w string. This S-matrix describes 
the scattering of the w part of an M\(v)w string with a y-particle. Alternately, fusing in u 
or v over the v roots making up a M\vw string gives 

M M M 

Y\ S 1 (u i -v)= j j S 1 (u-Vj) = SM+i(u-v)SM-i(u-v) = Y\_S 2 (ui-v), (D.3) 

which as we noted in the last equality, also coincides with fusion of the w-particle scattering 
matrix over the w roots making up an M\ (v)w string. By definition negative parity particles 
scatter with 

S (u-v) = S 1 (u + ik/g-v), (D.4) 

with positive parity particles, and regularly between themselves as we will come back to 
shortly. We would like to note that the fundamental scattering matrix for negative parity 
particles is almost inverse to the one for a k — 1 positive parity string, namely 

S (u-v)S k - 1 {u-v) = -l. (D.5) 

Next, to describe scattering of strings with one another this S-matrix is to be fused over 
the previously untouched argument. Note that there are two relevant cases for our consid- 
erations; either we fuse Sm over the v roots in a vw string, or we have to fuse Sm+\Sm-i 
over the w roots in a (v)w string. Of course either fusion gives the same result, as is 
immediately clear when interchanging the order of fusion, upon noting that the result is 
symmetric in the interchange of M and N. Let us define Smn as 

JV N M 

S M n(u-v)= jj S M (u-Vj)= [J |[ Si(ui-Vj), (D.6) 
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where the first product runs over the constituents of an N\(y)w string. This product can 
be rewritten as 

min(M,iV)-l 

S MN (U-V) = S M +n(u- V)S {M _ NI (U- V) Yl S \M-N\+2m( U ~ V ) ■ ( D - 7 ) 

m=l 

manifestly showing its symmetry under interchange of M and N. Negative parity particles 
scatter with positive parity strings with 

Som(u - v) = Si M {u + ik/g - v) = S^ +1 (u - v)Sj}_ l {u - v) , (D.8) 

while between themselves they scatter as positive parity particles meaning 

S o(u - v) = Su(u - v) . (D.9) 

Note that here we have 

S 0M {u-v)S k -i,M{u-v) = l. (D.10) 
Finally, note that these S-matrices are trivial when one of the indices is equal to k 

«Sfc = — 1, Smu = SkM = 1 • (D-11) 

This shows that at k these solutions to the discrete Laplace equation have a natural bound- 
ary. 

S yQ , S Qy and S^ M 

Fusing the scattering matrix of y ± particles with fundamental particles over a Q-particle 
bound state directly gives 



0/2 iu)- x ( v ) jx+{v) 



S™(«, v) = — J ^ , (D.13) 

_L_ _ x+{v) y x ( V ) 



where x are the parameters for a Q-particle bound state; x (v) = x(v ± iQ/g), and 
the subscript ± in S± denotes the sign of the imaginary part of the y-particle under 
consideration. Analogously we define the S-matrices for scattering of bound states with 
y-particles as 



cQv/.. „\ _ „Q/2 g (u)-x(v) x+{u) 

S - {u ' V) ~ q x+(«)-x(«)VFM' (D - 14) 
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As shown above in section B.2, the scattering matrix of a Q-particle bound state with an 
M\vw string of positive parity is given by 

x+(u) — x + (-y) x+(u) — x (v) x (u) 

where x =t (t') = x(t> ± iM/g), while for length one vw strings with negative parity we have 
simply 

qQOf., „\ = n Q ^ZM ~ x( - v + ^ k + jj/g) ~ g("+»(fc-l)/g) m 17l 

*W«,t;,_ g _ ^ + + 1)/g) x+{u) _ x~ (*u) ' 1 • ' J 

Again there is a special relation between these S-matrices, namely 

Sf£(u,v)S& k -Hu,v) = (-1) Q . (D.18) 

cQP 
°sl(2) 

In the present paper we will not derive the expression for the fused dressing factor, but 
assume structural analogy to the undeformed case and the effects of the deformation on 
the other kernels. Introducing the same split as in the undeformed case we write 

S S ^2)( X U X 2) = S 2 1 T,(x 1 ,x 2 y 2 , (D.19) 
where we have introduced the improved dressing factor E 

S(x 1 ,x 2 )= i ai " 2 q(x 1 ,x 2 ). (D.20) 
J- — + 

x l x 2 

We then assume these factors to fuse similarly to the undeformed case, to give 

-,QM _ ~-l v -2 
's[(2) ~~ °QM^QM 



S?i(2) — Sqm^qm ' (D.21) 



where now Tjqm should be a solution of the discrete Laplace equation in the bulk, and on 
the boundary at one with a source term analogous to the undeformed case [61]. 

D.2 Kernels and their properties 
K M and K M n 

As in the main text, we have 

and 

min(Af,JV)-l 

K M n{u) = Y^d^ l ° gSMN ^ = Km+N + K \M-N\ +2 ^ ^\M-N\+2j ■ (D.23) 

i=i 
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The kernels involving scattering of the negative parity particles can be defined in terms of 
the above via (D.5) and (D.10); note that we do not define them independently. For values 
of M greater than 2k we would have to take into account the 2k periodicity of the S-matrix 
and take M mod 2k in the above formulae, but we will not encounter this case. In order 
to simplify the TBA equations we would like to understand the integral identities satisfied 
by these kernels. 

For real rapidities the kernel Km is real and positive and has the following Fourier 
t r ctns fo i*ixi 

*„ ( „ ) = £^„ (U) = -E!^, (D , 4) 

where we defined the Fourier transform with an unconventional factor of g. The Fourier 
transform of the kernel Kmn is then 

K MN (u) = C °^ haJ (cosh (\M -N\- k)uj - cosh (M + N- k)u) - 6 M N ■ (D.25) 
sinn kco 

These kernels satisfy the following properties 

Kn{Sn,m - Inms) = s5 M ,i , , (D.26) 
Kml(Sl,n ~ Ilns) = sIm,n ■ , (D.27) 



where s(oj) = ^ — K— is the Fourier transform of s 



g f°° e —w*>u 



S{U) = L = I^sh^ ' (D - 28) 

and Imn = 5m,n-i + ^Af,iV+l is the incidence matrix, to be appropriately interpreted on 
the boundary as ijv,fc-i = ^Af,fe-2 34 - Since we will need it again, let us define 

(K + iy 1 = 1-1*8. (D.29) 

K yQ , K^, Kg M and 

In line with the conventions in the undeformed case, in the main text we defined 

Kg M (u,v) = ±,±lo g S^ v M (u,v) , (D.30) 

K%£(u,v) = -^LAi ogS ^( u>u)j (D.31) 

K$ y (u,v) = ^logS$ y (u,v), (D.32) 

1 d 

2iri du 



Kf(u, v) = f3— — log Sf{v, u) . (D.33) 



The scattering kernels for negative parity strings can be defined in terms of the above 
kernels via (D.18). We will mainly work with linear combinations of the last two kernels, 



Note that we have Kk = Kmh = 0. 
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namely 

K Qy (u, v) - K2 y (u, v) = K Qy (u, v) = K(u + iQ/g, v) - K(u - iQ/g, v) , (D.34) 

K Qy (u,v) + K2 y (u,v) =K Q (u,v), (D.35) 

K yQ (u,v) - Kf(u,v) = K Q (u,v) , (D.36) 

K yQ (u, v) + Kf(u, v) = K yQ (u, v) = K(u, v + iQ/g) - K(u, v - iQ/g) , (D.37) 

where 

w ) = ^ 3~ lo § H « ( D - 38 ) 

2m au x{u) — x(v) 

and Kq is defined in (D.22). These kernels satisfy the following properties under application 
of the operator (K + 

K xv N (3n,m - Inm * s) = 5q-i jM s + hi,iK Qy is, (D.39) 

k Zx($p,Q ~ Ipq * s ) = Sm+i,qs + 6 Qt iK M *s, (D.40) 

K yP (5p t Q- I PQ *s) = 5q, 1 (2K*s + s), (D.41) 

Kp{5p tQ -Ipq-ks) = 5 Qtl s, (D.42) 

where we note again that for M (not Q-particle) type indices ijvfc— l = while for 

Q type indices this incidence matrix is fine for Q = k — 1, but at Q = k we need slightly 
different identities. Namely 

K ^ix( 5 P,k ~ 2<5p iA: _i -ks) = 5 M ,k-is - KM,k-i * s , (D.43) 

K yP (5p, k - 25 Pifc _i * s) = . (D.44) 

The kernels entering in the above identities are 

1 d , x ( u ) 



K M = m -u b )--^^-^ y (D.45) 
and 

K M (u, v) = K(u + iM/g, v) + if(u - iM/g, v) . (DAG) 

Let us also define 

£ = log £±| . (D.47) 

In the above formulas x s is given by 



(e¥ sinh ^ - ^ sin 2 f ) 1 + 1 - ^^T7^ f 
;i + 5 2 sin 2 f)<? 2 sin 2 f 



sin 2 |' 
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K QP 

The main kernel K S U2) has the following structure 

= ~ K QP(u -v)- 2K% P (u,v) . (D.49) 

Proceeding by analogy to the undeformed case [4], we assume that the following identity 
[61] is valid 

K^ P -k (K + 1)^ = 6 ltR k$l8 , (D.50) 
where the kernel Kq,{u, v) vanishes for \v\ < ui,, and (K + is defined in (D.29). 

E The simplified TBA equation for Y k 

To derive the simplified TBA equation for we act with 8q & — 2Sq k—i s ° n the canonical 
TBA equations for Q-particles and make use of the special identities (D.43) and (D.44) 
and definitions (D.36) and (D.37) to get 

\ 2 

1 \ v(<*) 



log Y k = 2 log Y h . x * s + log ( 1 + Y^[ ]vw * 3 



k — 1 \vw 



log f 1 + -^-) * K MJk -i * a - log (l + YjM.) * * s (E.l) 

V Y M\vwJ V ' 

1 1 

log —^±2 * #fc-i * s + log (1 + Ip) * (K$ 2) - 2K™y 1 * s) . 



Note that the energy <?q is annihilated by the operator Sq^ — 2SQ,k-l8 because £k+i = 
£k-i and it satisfies the discrete Laplace equation by construction. Examining eq.(E.l) 
we immediately recognize a large part of the canonical TBA equation for Y),_\\ VVJ , albeit 
integrated with s and without the Yq term, for each a. Rewriting then gives 

logn=21ogr fc _i*s + log(l + -^— ) *s (E.2) 

V k-l\vw / 

+ log (1 + Y P ) * (K$ 2) - * s - 2K2~ l * s) . 

The kernel K s u 2 ) satisfies the discrete Laplace equation, meaning that if we introduce the 
would-be fused k + 1st kernel we should have 

K BI(2)- 2K S 1(2) * S - K S l(2)( K+1 )Q,k- K S l(2) * 8 + K ,t(2) * S -( K sI(2) - K sl(2) -0pk-l)*8 ■ 

In this way we reduce the equation for Yk to the following form 

lo g y*=21ogn_i*s + log('l + ^— *s-log(l + y fe _i)*s (E.3) 

\ k — 1 1 vw / 

+ log (1 + Y P ) * (K™! 1 - K^ 1 - 2K2- 1 ) * - • 
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Finally, the combination of the kernels in eq.(E.3) can be evaluated as follows. Let us 
consider the corresponding combination of S-matrices 

|pj^ («, v)) " = Sfa («, t, + f ) 5$, («, t, - f ) (S^-Hu, „)) ^ . 

Here we have used the fact that S^Ss is obtained by fusing the S-matrices S P ^ • Now we 
immediately realize that the first two terms are related by crossing. Using fusion in the 
first argument it is not difficult to show that 

^kUv + f)s^L(u,v-f)= (E.4) 



a pl (u, v + & V 1 (V « ~ f ) I"' = 9 2P /i P1 (u,v + f~ 2 

where a pi is the dressing factor and we have used the crossing equation describing scat- 
tering of a P-particle bound state with the fundamental particle. Here 

P Q , v x(u - iP/g) - x(v + iQ/g) 1 ~ x{u+iP/g)x(v+iQ/ 9 ) , . 

* (»•»' = ^-^i-x^-io/rti- ^,^., n s «-«<» - •> • 

as defined in appendix A. Thus, we are left to compute the quantity 

q- p h pl (u,v + ik/g)S™- 1 (u,v). (E.5) 

To proceed we will denote x(u ± iP/g) and x(y ± i(/c — l)/<?) as xf and x 2 • Note that 
x(v ±i/g + iA;/^) = x(u =F — l)/g) Tl = (x^) Tl . Working out the product (see also eq. 
(D.16)), immediately removing the shift on v in h as just indicated, we get 



q y h^ L (u,v + ik/g)S P v (u,v) 
x7-4 1 



1 ^2 X l _ X 2 X l _ X 2 X l 

1 x 2 1 £+3,+ x l x 2 x l x 2 x l 



1 ^2 



P-l fc-2 

x Si-p +2 j{u,v + ik/g) 5 P+ i_ fc+2 m(w, u ) 

j'=l m=l 



1 P-l fc-2 



- Xl i 2 X \ X + ]J Si- P+2 j(u,v + ik/g) Y[ Sp+i-k+2 

J. _i_ _i_ X i 3^ i 



cc^a;^ 1 2 i=l m=l 

P-l fc-2 

g -CP-i)5 p+1 _ fc ( U}U ) S^p^jfav + ik/g) Sp+ X - k+2m (u,v) 

j=l m=l 

-1)V (P_1) . (E.6) 



In the third equality we used the general identity 

1 



V + _ i = q- {P+Q) S P - Q , (E.7) 
xtxf X i X 3 
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where P and Q refer to the bound state numbers of i and j respectively. In the fourth we 
used 

M-l 

Y[ 5i-P+2j = Sp-i , (E.8) 

3=1 

S P (u±ik/g) = -Sp_ k (u) (E.9) 



and 



k 



n^=(-i) p+fc+i . (e.io) 

3=1 

Thus, we have proved that the combination of the S-matrices we are interested in is constant 
and therefore the corresponding combination of the kernels vanishes. With this result eq. 
(E.2) simplifies further, now properly deserving the name, giving 

lo g y fe = 21ogY fc _ 1 * S -log(l + y fc _ 1 )* S + log(\ + F ^ N ) *s. (E.ll) 

\ k — 1 1 vw J 

Note that the contribution of Yj,_\\ vw has doubled up here, exactly as it did in the equations 
for vw strings, though admittedly in a less obvious fashion. 

F Representations of su 9 (2|2) at roots of unity 

Here we discuss some elements of representation theory of the centrally extended su g (2|2) = 
Uq(su(2\2)) at even roots of unity. We follow the construction of [71]. 

Quantum oscillators 

Let us first introduce quantum oscillators. We will use those to build our representations. 

The g-oscillator (q-Heisenberg-Weyl algebra) U q (t)4,) is the associative unital algebra 
consisting of the generators {a^, a, w, w -1 } that satisfy the following relations 

w a^ = q a'w, qwa = aw, (F-l) 

ww~ l = w~ 1 w = l ) a a^ — q a^a = w" 1 . 



-\ah 

such, we will set it to zero in the remainder. The defining relations then imply that 



From the defining relations one can see that the element w 1 (a^a — w q _™-i ) is central. As 



ata = "~\ aa^ = qW - q ^\ (F.2) 

q-q~ l q-q- 1 

Since we are interested in representations of supergroups, we will also need to consider the 
fermionic version of the g-oscillator. The above notion is extended to include fermionic op- 
erators by adjusting the defining relations in the following way (we keep the same notation 
for bosonic and fermionic a, at for now) 

wa^ = qa^w, qwa = aw, (F.3) 

ww^ 1 = w~ 1 w = 1, a a^ + qa^a = w. 
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resulting in the following identities 



In this case, the central element is w (ata - ts ^=r). A gain we set this element to zero, 



*U = W - W ~" ^ = qW ~ l -\ lw . (FA) 

q-q q-q 

Of course in the fermionic case the operators a, a^ square to zero. From (F.4) it then follows 
that this only is consistent if w 2 = l,q 2 . Below we will identify w = q N , where N = 0,1 is 
the number of fermions making it indeed compatible. 

Fock space 

Let us first build the Fock representation of Uq(f)i). For this purpose consider a vacuum 
state |0) such that 

a|0)=0, H0) = |0), (F.5) 

then the Fock vector space T generated by the states of the form 

\n) = (at)«|0) , (F.6) 

is an irreducible module of U q (t)4). Let us first consider the bosonic (/-oscillators. With 
the help of the defining relations (F.l) and (F.2) one finds that the action of the oscillator 
algebra generators on this module is 

a'|n) = |n + 1), a|n) = [n] q \n — 1), w\n) = q n \n). (F.7) 

This makes it natural to identify w = q N , where N is understood as a number operator. 
Analogously, fermionic generators are found to act as 



a' Ira) 



\n + l), a|n) = [2 — n] q \n — 1), w\n) = q n \n). (Fi 



However, due to the fermionic nature, n can only take the values and 1 and thus in (F.8) 
we can identify [n — 2] q with [n] q . 

Fock space at roots of unity 

Let q = e~T* 35 ; then it is easy to see that a 2fc , (a^) 2k , w 2k , w~ 2k are central elements. Of 
course, for fermionic oscillators, this is trivial and holds always. Because of this, the Fock 
space is reducible and has a finite number of weight spaces. Indeed, since 

a|rj) = [n] q \n — 1) = 0, w\n) = q n \n) (F.9) 

we see that only the weights {1, . . . , q 2 ^ 1 } are present. In other words, we can restrict our 
states in the Fock space to n = 0, . . . , k — 1 only and doing so gives a irreducible module 
again. 



35 All the results are straightforwardly generalized to the case where q = e k , with k an odd integer. 
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This allows us to introduce the class of so-called semi-cyclic representations of the 
quantum oscillator algebra. Let (A,/x) G (C x {0}) U ({0} x C), then we define 



a|n) = [n] q \n- 1), for n = 1, . . . k - 1 a|0)=/z|fc-l) (F.10) 

a t| n ) = |ra + l), for n = 0, ...k-2 a f |fc - 1) = i/|0) (F.ll) 

w|n) = g n |n), for n = 0, . . . k - 1. (F.12) 

The idea behind this is easy; we identify \k) ~ |0). Putting fi = v = simply returns us to 
the normal restricted Fock space. 

Representations of U q (5l(2)) 

We will now construct representations of U q {s\{2)). Consider two copies of bosonic q- 
oscillators a^aLtu, = q Ni which mutually commute. Then the Fock space is naturally 
spanned by vectors of the form 

\m,n) = (a\r(4r\0). (F.13) 

Define the subspace Tm = span{ \m, M — m) \ m = 0, . . . , M }. We will refer to M as 
the total number of particles and dimJ-M = M + 1. It is easy to see that under the 
identification 

£ = a|ai, F = a\z 2 , H = N 2 - N u (F.14) 

we obtain a presentation of C/ g (sl(2)). Then depending on the representation of the q- 
oscillators on the Fock space (normal or semi-cyclic) this will define an explicit represen- 
tation. 

Irreducible representations 

Generic q. Let q be a complex number that is not a root of unity. In this case the 
q-oscillators have the standard representation on the Fock space (F.7). Moreover, the sub- 
space Tm = span{ \m, M — m) \ m = 0, . . . , M } is an irreducible £Y,j(s [(2) ^representation 
of dimension M + 1. Explicitly we have 

E\n, M — n) = [n] q \n - 1,M - n + 1), F\n, M — n) = [M - n] q \n + 1,M - n - 1), 
H\n, M - n) = (M - 2n)\n, M — n). (F.15) 

It is easy to check that indeed satisfies all the defining relations of U q (si(2)). 

Root of unity q = e~k . Now we have to use the semi-cyclic presentation of the algebra. 
Let //i, v\ and /i2, ^2 be the parameters describing the representation of the first and second 
set of oscillators respectively. Notice that the only when the total number of particles is 
k — 1 there is a possibility of seeing the semi-cyclicity. Thus, when restricting to states 
with total M < k — 1 number of particles, representations coincide trivially with those for 
generic q given in (F.15). 
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However, when M = k — 1 we find a one-parameter family of representations. This is 
obtained by taking fJ,x = u 2 = 0. Then apart from the ordinary relations for n = 1, . . . , k — 1 
(F.15) we find 

E\0,k- 1) = ^iM 2 |&- 1,0). (F.16) 

Alternatively, we can set 1^1 = ^2 = and get a special relation for F acting on |0, k — 1). 

The special point where fa = Vi = is most natural in the sense that the represen- 
tations simply correspond to the restriction of the normal q-oscillator representations like 
for M < k — 1. All these representations are irreducible. 

Reducible representations at roots of unity 

By construction, any representation based on the normal oscillator representation on Tm 
with M > k — 1 for q = e~fc" is reducible. At M = k, the representation contains two singlets 
|0, k), \k, 0). In general, the representations will be indecomposable. More precisely, let 
M = nk + m with < m < k then the representation is reducible and contains 



0^ m CJ M . (F.17) 



i=0 

Only when m = k—1, the representation is decomposable. The above can be readily proven 
by using that the copies of T m are generated by \ak, M — ak), \ak + 1,M — 1 — ak), . . . 
where a = 0, . . . , n. 

Representations of centrally extended su g (2|2). 

We will now construct the bound state representation for centrally extended su,j(2|2) in 
the g-oscillator language. We need to consider two copies of C/q(sl(2)), a bosonic and a 
fermionic one. Thus we need four sets of g-oscillators a^, a\,Wi = q Ni , where the index 
i = 1, 2 denotes bosonic oscillators and i = 3, 4 - fermionic ones. Using these we write 

£i = ajai, F^al^, H 1 =N 2 -N 1 , (F.18) 

+ + + + N-i + — No — Na 
E 2 = a a[a 2 + b a[a 3 F 2 = c a^i + d aJ,a 4 , H 2 = -C+— -, (F.19) 

£ 3 = aja 4 , i^alas, H 3 = N 4 - N 3 , (F.20) 

where C is central. Let V = q c . We first consider the case where q is not a root of unity. 
It is then straightforward to check that this set of generators forms a representation of 
su g (2|2) on the Fock space when restricting to the subspace of total particle number M 
upon setting 

. . [C ~ , 1 ~ U 2 V 2 , 9 V- 2 ~ U- 2 rp9n 

ad = — - — ^ — , be = — - — ^ — , ab = qa — - — , , cd = — : . r.21 

[M] q [M] q ' y [M] q a [M] q V 1 

The dimension of this representation is AM. 
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Let us spell out the representation more explicitly. The bound state representation is 
defined on vectors 



Kn,M>H4n4)Vi) fe (4)'|0>, (F.22) 

where the indices 1, 2 denote bosonic oscillators and 3, 4 - fermionic and the total number 
of excitations k + l + m + n = M is the bound state number. 

The triples corresponding to the bosonic and fermionic sl 9 (2) in this representation 
are given by 

Hi\m, n, k, I) = {I — k)\m, n, k, I), H^m, n, k, I) = (n — m)\m, n, k, I), 

Ei\m, n, k, I) = [k] q \m, n, k — 1, 1 + 1), E^\m, n, k, I) = \m + l,n — 1, k, I), 
Fi\m, n, k, I) = [l] q \m, n, k + 1, 1 — 1), F^\m, n, k, I) = \m — 1, n + 1, k, I). (F.23) 

The supercharges act on basis states as 

H2\m,n,k,l) = — |(7 1 \m,n, k,l), 

Ei\m, n, k,l) = a (— l) m [/] 9 \m, n + 1, k, I — 1) + b \m — 1, n, k + 1, 1), 

F 2 \m, n, k,l) = c [k] q \m + 1, n, k - 1, 1) + d (-l) m |m, n—l,k,l + 1). (F.24) 

For generic q these representations are irreducible. 

Root of unity. At the root of unity, we can once again use the semi-cyclic representation 
of the quantum oscillators. This will again yield a special representation at M = k — 1, 
which upon setting \±i = v i = reduces to the one coming from the normal representation 
on the Fock space. 

Reducible representations at roots of unity. In the case where [n = v i = one 
finds reducible representations for M > k. Since this representation contains a bosonic 
su(2), the discussion of the reducibility of those representations also applies here. In fact, 
we immediately see that we get indecomposable representations again. Let Vm be the 
M-particle bound state representation and let M = nk + m with < m < k. Then, it is 
not difficult to see that we again find 



®K«C%. (F.25) 



i=0 



The central charges obviously are the same in all components and consequently each of the 
sub representations share the same central charges U, V . 
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